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yjQ ' Abstract. In this paper we have analyzed 25 Galactic O and early B-stars by means of H and A* band spectroscopy, with the 

primary goal to investigate to what extent a lone near-IR spectroscopy is able to recover stellar and wind parameters derived 
\Q ■ in the optical. Most of the spectra have been taken with subaru-ircs, at an intermediate resolution of 12,000, and with a very 

high S/N, mostly on the order of 200 or better. In order to synthesize the strategic H/He lines, we have used our recent, line- 
blanketed version of fastwind iPuls et aljEool . In total, seven lines have been investigated, where for two stars we could 

in . 

make additional use of the Hei2.05 singlet which has been observed with irtf-cshell. Apart from Br y and Hen2.18, the other 
lines are predominately formed in the stellar photosphere, and thus remain fairly uncontaminated from more complex physical 
processes, particularly clumping. 
First we investigated the predicted behaviour of the strategic lines. In contradiction to what one expects from the optical in the 
O-star regime, almost all photospheric H/Hei/Hen H/K band lines become stronger if the gravity decreases. Concerning H and 
Hen, this finding is related to the behaviour of Stark broadening as a function of electron density, which in the line cores is 
different for members of lower (optical) and higher (IR) series. Regarding Hei, the predicted behaviour is due to some subtle 
, NLTE effects resulting in a stronger overpopulation of the lower level when the gravity decreases. 

We have compared our calculations with results from the alternative NLTE model atmosphere code cmfgen iHillier & MiiieJ 
1998). In most cases, we found reasonable or nearly perfect agreement. Only the Hei2.05 singlet for mid O-types suffers from 
, some discrepancy, analogous with findings for the optical Hei singlets. 

For most of our objects, we obtained good fits, except for the line cores of Br y in early O-stars with significant mass-loss. 
Whereas the observations show Br r mostly as rather symmetric emission lines, the models predict a P Cygni type profile with 
strong absorption. This discrepancy (which also appears in lines synthesized by cmfgen) might be an indirect effect of clumping. 
After having derived the stellar and wind parameters from the IR, we have compared them to results from previous optical 
analyses. Overall, the IR results coincide in most cases with the optical ones within the typical errors usually quoted for the 
corresponding parameters, i.e, an uncertainty in T e g of 5%, in logg of 0.1 dex and in M of 0.2 dex, with lower errors at higher 
wind densities. Outliers above the 1-cr level where found in four cases with respect to logg and in two cases for M. 

Key words. Infrared: stars - line: formation - stars: atmospheres - stars: early type - stars: fundamental parameters - stars: 
winds, outflows 



1. Introduction 

Although rare by number, massive stars dominate the life cycle 
of gas and dust in star forming regions. They are responsible for 
the chemical enrichment of the ISM, which in turn has a sig- 
nificant impact on the chemical evolution of the parent galaxy. 
The main reason for this is that due to their large masses, each 



physical stage evolves on much shorter timescales and more 
violently, compared to low-mass stars, which provides a very 
efficient recycling of elements. Moreover, the large amount of 
momentum and energy input of these objects into the ISM con- 
trols the dynamical evolution of the ISM a nd, in turn, the evo- 
lution of the parent galaxy (e.g..lLeitherer & Heckman II 19951 
ISilich & Tenorio-TaglJ200ltlOevl2003l) .~ 
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Presently, high mass star formation is still poorly under- 
stood. By their nature, all star forming regions are found buried 
in molecular gas and dust, allowing little or no light to escape 
at optical wavelengths. But, unlike low mass stars which may 
eventually become optically visible while still contracting to 
the main sequence, the short contraction timescale of high mass 
stars keeps them deeply embedded throughout the entire for- 
mation process. Compared to the optical, the dust and gas sur- 
rounding young massive stars become more transparent in the 
infrared (IR) regime. Observations at such wavelengths reveal 
the hot stellar content of these dust-enshrouded environments 
like young HII regions in dense molecular clouds, the Galactic 
centre or massive clusters. Recent examples which illuminate 
the advantages of near- and mid-IR wavelengths for observ- 
ing massive proto-stars and star formation triggered by mas- 
sive st ars h ave been given bvlBlum et al.l d2004 . 1whitnev et alJ 
d2004 and lClark & Porter! d2004 ~ 

Following the substantial progress in ground-based IR in- 
strumentation in the past decade, IR spectroscopy has become 
a powerful diagnostics for the investigation of hot stars and 
the stellar winds surrounding them. The first systematic obser- 
vational studies of OB stars in the H a n d K band have been 
perf ormed bv. e.g.. lHanson et all dl996l) . iMorris et ail dl996l) 
and iFullerton & NaiarroT(ll998i) ~providing an important basis 
for quantitative spectral analysis of early type stars. With the 
use of satellites (e.g., the Infrared Space Observatory (ISO) in 
1995 and the Spitzer Space Telescope in 2003) a larger spectral 
window became accessible, completing the IR regime already 
observed from the ground. 

Modeling of the near-infrared, on the other hand, has been 
performed mostly for early-type sta rs with dense winds, i.e., 
for Wolf-Rayet Sta rs dHillierlll982l). Galactic centre objects 
dNaiarro et alJl994 . Of/WN stars dCrowther et al J 1 995ll 1 9981) 
and Luminous Blue Variables dNaiarro et al.l 119971 I1998I) . 
Combined ground-based optical and near-IR, ISO and Spitzer 
mi d-IR spectra for LBV s and Wolf-Ra yets ha ve b een mod eled 
by Naiarroet alJll 997hUDessart et al JfcOOdl andlMorris et all 
(2000, 2004), in a similar spirit as outlined below, namely to 
compare optical parameters with those obtained over the wider 
spectral range, partly including also the UV. Note that all these 
investigations have been performed by means of the model at- 

J — I 73 

mosphere code cmfgen (Hillier & Miller 1998). 

For objects with thinner winds (which are of particular 
interest when aiming at the youngest objects emerging from 
Ultra-Compact Hn (UCHII) region s), no results ar e avai lable 
so far, except from a pilot study bv iLenorzer et alJ d2004 . In 
this study, various synthetic H/He IR-profiles, located in the J 
to L band, are presented for a comprehensive grid of O-type 
stars (from dwarfs to supergiants), and their diagnostic poten- 
tial and value is discussed. 

The reader may note that most of the available datasets 
of IR-spectra have been observed at relative ly low resolu- 
tion (typically, at R « 2,000, though Fulle rton & Nai arro 
(1998)' present a few spectra with R « 10,000), which com- 
promises a precise spectroscopic analysis, since many deci- 
sive spectral features rem ain unresolved. Meanwhile, how- 
ever, lHanson et al.1 d2005l) have re-observed a large sample 
of Galactic O-type "standards" with much higher resolution, 



typically at R » 12,000. The objects where chosen in such 
a way that they both largely overlap wit h stars which have 
been analyzed befor e in the optical (e.g.. iHerrero et ali r2002: 
iRepolust et alJEoO^) . and cover a wide range in spectral type 
and luminosity class. Therefore, the present paper has the fol- 
lowing objectives: 

• We carry out a spectral analysis for this sample in the near 
infrared regime and compare it with results already ob- 
tained in the optical. This will allow us to check the extent 
to which the data derived from the IR is consistent with 
results obtained from alternative studies in different wave- 
length bands. As an ultimate goal, we plan to use solely 
the infrared regime to provide accurate constraints to the 
characteristics of stars which can only be observed at these 
wavelengths. 



We 


test our model 


atmosphere 


code 


FASTWIND 


(Santolava-Rev et al. 


1997; 


ITT — — n 

Herrero et al. 


2002; Puis et al. 


2005 


) for OB stars in the near infrared, hence extending its 



usage to these wavelength ranges. 
• We give special attention to those lines which are located 
in the H and K band, i.e., which can be accessed by 
ground-based instrumentation alone. Note that these lines 
are mainly formed close to the photosphere, i.e., remain un- 
contaminated by additional effects such as clumping and 
X-rays and, thus, should provide rather robust estimates for 
effective temperatures and gravities. 

The remainder of this paper is organized as follows. In 
Sect. |2]we briefly describe the observations and the lines used 
in our analysis. In Sect. [3] we summarize our model calcula- 
tions and comment on our treatment of line-broadening for the 
hydrogen lines. Sect. 0] outlines some theoretical predictions 
concerning the behaviour of strategic li nes, and Sect. Ecom- 
pares our results with those obtained bv lLenorzer et al.l (12004) 
by means of the alternative wind-code cmfgen (see above). In 
Sect. [6] we discuss the analysis of the individual objects of our 
sample, and in Sect. we consider the consequences related to 
missing knowledge of stellar radius and terminal wind velocity. 
Sect.|8]compares our IR results with those from the correspond- 
ing optical data. In Sect.|9] finally, we present our summary and 
conclusions. 

2. Observations, targets and strategic lines 

For ou r analysis we use a subset of stars given bv lHanson et alJ 
(2005). Detailed information about the observation dates, res- 
olution, spectrometers and data reduction can be found there. 
We selected the spectra from that sample which where obtained 
with the Infrared Camera and Spectrograph (ircs) mounted at 
the Cassegrain focus of the 8.2m Subaru Telescope at Mauna 
Kea, Ha waii. This totaled in 2 9 stars out of the 37 targets col- 
lected bv lHanson et al.1 d2005l) . 

The targets had been selected i) to adequately cover the 
complete OB star range down to B2/B3 at all luminosity 
classes, and ii) that most of t hem have already b een analyzed 
in the optical (for details, see lHanson etal1 l2005). According 
to the purpose of our analysis, we have exclusively used the 
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Table 1. Sample stars and observing data in the H and K band. 
In sub-samples I to III we have grouped those objects which 
have been previousl y analyzed in the optical (sub-sample I: 
iReoolust et al.ll2004 sub-sample II: [Herrero et alJl2000ll2002t 
sub- sample III: lKudritzkietaDll999l) . Subsample IV com- 
prises those objects covered by various authors or not analyzed 
at all. 
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only K band available. 

Additional irtf-cshell spectra covering Hei2.05 available. 



data from the Subaru Telescope and not the VLT data (com- 
prising the remaining 8 objects), since we did not possess com- 
plementary optical spectra for the latter dataset. In the follow- 
ing, we will define four different sub-samples denoted by I 
to IV in order to distinguish between objects analyzed in the 
optical by dif ferent authors . Samp le I comprises those stars 
discussed by Repolust et al. (2004), s ample II c orresponds to 
objects analyzed bv iHerrero et al.l yOOO, 200 21) sample III 
(B-supergiants) has been analyzed bv iKudritzki et all i l 19991 
only with respect to wind-parameters), and sample IV con- 
sists of the few remaining objects considered by various au- 
thors or n ot at all. In p a rticula r, HD 46150 has been inves- 
tigated by IHerrero et alJ dl992l plane-parallel, unblanketed 

1 Note that the first of the two investigations has been performed by 
unblanketed models. 



models) and r Sco (HD 149438) bv lKilian et all dl99ll plane- 
parallel NLT E analysis w ith underlying Kurucz models) and by 
IPrzvbilla & Butlen(l2004l) with respect to optical and IR hydro- 
gen lines. Table^indicates to which individual sub-sample the 
various objects belong. 

The Subaru/iRcs H band and K band spectral resolution is 
R « 12000. The typical signal-to-noise ratios obtained with 
these spectra were S/N » 200-300, with areas as high as S/N 
~ 500, and as low as S/N * 100, depending on the telluric con- 
tamination. The spectra were obtained over two separate runs, 
the first in November 2001 and the second in July 2002. Due 
to poor weather condition, the telluric corrections for some of 
the spectra proved to be difficult. This can be seen in the H 
band spectra of HD 217086, HD 149757, HD 66811, HD5689 
and HD 15629. Furthermore, there were no H band spectra of 
HD 15558 and r Sco available, weakening the significance of 
their analyses. The reduction of the data was performed using 
iraf routines and Perl idl including standard procedures such as 
bias subtraction, flat field division, spectrum extraction, wave- 
length calibration and continuum rectification. TableHsumma- 
rizes all observational runs obtained with ircs. In the following, 
all wavelengths of NIR lines are given in microns (//m). 

The data for the Hei/i2.05 line, which had not been ob- 
served by subaru were taken at the Infrared Telescope Facility 
(irtf) in Mar ch, June and July o f 2003. The cshell echelle 
spectrograph iGreene et alJll993l) was used with a slit of 1.0 
arcseconds. The instrumental spectral resolving power as mea- 
sured by a Gaussian fit to the OH night sky emission lines was 
4.0 pixels FWHM, or 12.1 km s _I , corresponding to a resolu- 
tion of 24,000. The spectra were reduced using iraf routines 
and the subsequent analysis was done using routines written 
in Perl idl. For all spectra, dark frames and flat field frames 
were averaged together to form a master dark and flat frame. 
Unfortunately, Hei2.05 lies within a region where the telluric 
absorption is extremely large , degrading the signal significantly 
dKenworthv & Hansonl2004T) . After the reduction, it turned out 
that most of our spectra did not posses sufficient quality (only 
moderate S/N), and we could use only the spectra obtained 
for two of the stars (HD 190864 and t Sco) for our analysis. 
Nevertheless, in all cases we have included the synthesized line 
for the sake of completeness. 

The spectral classification of sample I is the on e adopted 
bvlHerrero et al.l {l992), based mostly on the work bv lWalbornl 
dl972l Il973l) , the unpublished catalogue of OB stars by 
C. Garmany and by Mathvs (1989). As for sam ples II to IV, the 
spectral classification used by Hanson et al ] J2005I) has been 
retained. The classifications were based mostly on Walborn 
classifications, except for th e Cyg OB2 stars, which relied on 
Mas sev & Thompson! II 1991). 

In total the sample consists of 29 Galactic O and early B type 
stars as listed in Tablefflranging from 03 to B3 and covering lu- 
minosity class la/lab, Ib/II, III, and V objects, where 4 stars (of 
the latest spectral types) have been discarded later in the study. 
The strategic lines used in our analysis are (all wavelengths in 
air) 

• //band 

- Hi/U.68(n = 4-»ll,Brll), 
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- HiA1.74(n = 4->10,BrlO), 

- Hei ,11.70 (3p 3 P° - Ad 3 D, triplet), 

- HenA1.69(n = 7 ->12). 
• T^band 

- Hi ,12.166 (n = 4 ^7, Br y ), 

- Hei ,12.058 (2s [ S - 2p 'P°, singlet), where available, 

- Hei ,12.1 1 (comprising the Hei triplet XI. 1 120 (3p 3 P° - 
As 3 S) and the Hei singlet ,12.1132 (3p 'P - As [ S)), 

- Hen ,12.188 (n = 7 -»10). 

Note that Br r overlaps with the Hei triplet ,12.1607 (Ad 3 D - I f 
3 F°), the Hei singlet XL, 1617 (Ad 'D - If l F°) and Hen ,12.1647 
(n = 8 — >14). Whereas the singlet is not included in our formal 
solution, the Hei triplet, in particular, has been used to check 
the consistency of our results. Note that the influence of the 
Hen lines overlapping with BrlO and Brl 1 is marginal. 

3. Model calculations 

The calculations presented in this paper have been performed 
by means of our present version of fastwind, as described by 
Puis et al. (2005). In addition to the features summarized in 
Repolust et alJ J2004I) . this code meanwhile allows for the cal- 
culation of a consistent 2 temperature, utilizing a flux-correction 
method in the lower atmosphere and the thermal balance of 
electrons in the ou ter one. As has been discussed, e.g., by 
iKubat et alJ l[l998), the latter method is advantageous com- 
pared to exploiting the condition of radiative equilibrium in 
those regions where the radiation field becomes almost inde- 
pendent on T e . Particularly for IR-spectroscopy, such a con- 
sistent T-stratification is of importance, since the IR is formed 
above the stellar photosphere in most cases and depends (some- 
times critically) on the run of T e . We have convinced ourselves 
that our previous results concerning optical lines remain (al- 
mo st) unaffec t ed by this modification. 

IPuls et al1ll2005h present a thorough comp arison with mod- 
els fr om alternativ e "wind-codes" (wM-ba sic, IPauldrach et alJ 
2001 and cmfgen, iHillier & Miilelll998h . Some differences 
were seen in the On continuum at and below 350 A (fastwind 
predicts a higher degree of line-blocking in this region), which 
might have some influence on the helium ionization balance, 
due to a different illumination of the Hen resonance lines. Also, 
cmfgen predicted weaker optical Hei singlets in the tempera- 
ture range between 36,000 to 41,000 K for dwarfs and between 
31,000 to 35,000 K for supergiants. Otherwise, the comparison 
resulted in very good agreement. 

3.1. Atomic data and line broadening 

In order to obtain reliable results in the IR, our present H and 
Hen models consist of 20 levels each, and Hei includes lev- 
els until n — 10, where levels with n = 8,9, 10 have been 
packed. Fu rther informatio n concerning cross-sections etc. can 
be found in l.TokufhvN2002h . 

The hydrogen bound-bound collision strengths require 
some special remarks. The atomic data on radiative line pro- 

2 Note, however, that non-radiative heating processes might be of 
importance. 



cesses in Hi are very accurate because they can be obtained 
analytically due to the two-body nature of the hydrogen atom. 
However, for excitation/de-excitation processes, these involve 
a colliding particle, making the situation much more complex. 
In most cases only approximation formulae are available. 

Note that the "choice" of the collisional data is an es- 
pecially important factor for the line formation in the IR. 
Although the effect of different collisional data will not be 
apparent for the ground state, higher levels display a signif- 
icant sensitivity, reaching its maximum fo r levels with inter- 
media te n at line formation depth. Recently, Przvbill a & ButleJ 
(2004) have emphasized the differences in the collisional cross 
section from approximation formulae and ab initio computa- 
tions for transitions up to n - 7. Particularl y, the fr eque ntly 
used approximations bv lMihalas et alJ i 19751) and bv Ijohnsonl 
dl972l) show a different behaviour and fail to simultaneously 
reproduce the optical and IR spectra over a wide parameter 
range. However, the collisional data provided by Przybilla & 
Butl er (in combinat i on wi th t he approximation for mulae by 
IPercival & Richards! Jl978h and lMihalas eTall dl975h are able 
to reproduce the observed line profiles in those cases which 
have been checked. Note, however, that these checks did not 
cover O-type supergiants, cf. Sect. 18. II 

The standard implementation of the corresponding cross 
sections i n fastwmd, on th e other hand, is based on data pre- 
sented by Giovan ardi et al.l (Tl987). Although affected by simi- 
lar problems as described above, the differences to the ab ini- 
tio calculations are smaller but still worrisome. As detailed 
later on, a comparison of simulations using both data-sets al- 
ternatively revealed that for our O-star sample we find bet- 
ter agreement with corresponding optical results if our stan- 
dard implementation is used. Consequently, all calculations de- 
scribed in the following are based on these data, whereas fur- 
th er comments concerning th e effect of incorporating the data 
bv lPrzvbilla & Butlerl (Eo04) are given in Sect.lO 

Since we are concentrating on those lines which are formed 
close to the photosphere, line-broadening is particularly impor- 
tant (and leads to a number of interesting effects, shown below). 
Unfortunately, calculations as "exact" as for optical lines do not 
yet exist for their IR counterparts, leaving us to use reasonable 
approximations. 

Actually, ll.emkel dl997l) has published extended 
Stark broadening tables (based on the approach by 
IVidal. Cooper & Smith! Il973l "VCS") for the hydrogen 
Lyman to Brackett series. In a first step, we have used his 
dataset for the calculation of the Brackett lines. However, we 
immediately realized that at least Brl 1 must be erroneous by 
virtue of a comparison with observed mid-resolution B-type 
NIR spectra which revealed no proble ms if app roximate 
broadening functions are used (see lHanson et al J20031 Fig. 4). 
After a careful investigation by K. Butler (priv. comm.), it 
turned out that not only Brll but also other transitions, i.e., 
predominantly members of the higher series, are affected by a 
number of (numerical) problems in the code used by Lemke. 

Thus, Stark-broadening of hydrogen needs to_be approx- 
imated as wel l. We follow the met hod by Grieml d 1967b 
as outlined in lAuer & M ihalas (1972, Appendix), based on 
a corrected asymptotic Holtsmark formula. Due to com- 
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Fig. 1. Comparison of strategic NIR lines for two atmospheric 
models at T^— 40,000 K and different gravities, logg=3.7 
(solid) and log g=4.5 (dashed), respectively. Both models have 
a negligible wind, with log Q = -14. The number in the up- 
per left corner gives the equivalent width (in A) of the low- 
gravity model, where, in agreement with previous papers, nega- 
tive numbers indicate net-absorption. All profiles are displayed 
on the same horizontal scale (of width 0.02yum), and the profiles 
have been rotationally convolved with V r sin /= 80 km s _1 . 



parisons with VCS calcu lations for optical transitions from 
ISchoning & Butle11 Jl989l) . which are used by fastwind any- 
way, we have convinced ourselves that the Griem approxima- 
tion recovers the more exact VCS results with very high pre- 
cision, if the upper and lower level of the transition lie not too 
closely together (e.g., H a is badly approximated, whereas for 
H r no differences are visible). The results obtained by using 
either the (erroneous) data by Lemke or the Griem approxi- 
mation are given by means of a detailed comparison later on, 
cf. Fig.|6] Griem-broadening is also applied to Hen (1.69, 2.18 
/im), whereas for Hei (1.70, 2.05, 2.1 1 jjm) we have used Voigt 
functions only, with damping parameters accounting for natu- 
ral and collisional broadening. The comparison to observations 
suggests that this approximation describes reality sufficiently 
well. 



4. Predicted behaviour of strategic lines 

Before we describe the results of our analysis, we will in- 
vestigate the behaviour of our synthesized lines in some de- 
tail, particularly because their dependence on gravity seems to 
be somewhat strange, at least if one extrapolates the knowl- 
edge accumulated in the optical. Although a related investiga- 



tion has already been performed bv lLenorzer et alJ J2004). they 
have only discussed the behaviour of the equivalent widths. 
Moreover, their model grid is rather restricted and does not al- 
low the investigation of changes in synthesized profiles if only 
one atmospheric parameter is altered. On the other hand, we 
have calculated a rather large grid of models in the parameter 
range 20, 000 < < 50, 000 with a typical variation in log g 
over two dex, and wind stren gths varying from negligible to 
very large (cf. Puis et al. 2005), allowing us to inspect this kind 
of reaction in more detail. In the next section we will, of course, 
compa re our results also to those obtained by iLenorzer et alJ 
d20Q4t) . 

As a prototypical example, in Fig ^ we compare the strate- 
gic H/He NIR lines for a model at T eS = 40000 K, with logg= 
3.7 (solid) and 4.5 (dashed). Both models have a vanish- 
ing wind density, corresponding to log Q = -14, where Q 
is a suitable measure to compare the in fluence of differe nt 
wind strengths in recombination lines (see lPuls et alJ Jl9 96)). 
Throughout this paper, we have defined 



M[M Q /yr] 



((RJR^v^kms- 1 ]) 1 



(1) 



Most interestingly, almost all NIR features become deeper 
and their equivalent width increases if the gravity decreases. 
In contrast to the Balmer lines, the cores of BrlO, Brll (and 
of Henl.69/2.18) are significantly anti-correlated with grav- 
ity. This behaviour is completely opposite to what one expects 
from the optical. Only the far wings of the hydrogen lines bear 
resemblance to the optical, which become shallower when the 
gravity decreases. 

Although the reaction of Hei on \ogg is only moderate, 
at lower temperatures (with more Hei present) we observe 
the same trend, i.e., the equivalent width (e.w.) increases with 
decreasing gravity, as shown in the e.w. iso-contour plots in 
Fig. |3 For comparison, this plot also shows the extremely 
"well-behaved" Hei447 1 line, which decreases in strength with 
decreasing gravity in all regions of the r e ff-log g -plane. 

Before we will further discuss the origin of this peculiar 
behaviour of NIR-lines, let us point out that these trends do 
not depend on specific details of the atmospheric model, par- 
ticularly not on the presence or absence of a temperature in- 
version in the upper photospheric layers. The same relations 
(not quantitatively, but qualitatively) were also found in mod- 
els with a monotonically decreasing temperature structure in 
the inner part (log tr oss > -2) and a constant minimum tem- 
perature in the outer wind. 

4.1. Hydrogen and Heu lines: Influence of Stark 
broadening 

The peculiar behaviour of the line cores of the hydrogen 
Brackett lines and Henl.69/2.18 can be understood from the 
reaction of the core of the corresponding Stark-profiles as a 
function of electron density. Fig. [3] shows the Stark-profiles 
for H r and BrlO as a function of frequency displacement from 
the line centre in units of thermal Doppler-width, calculated 
in the Griem approximation. Both profiles have been calcu- 
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Fig. 2. Iso-contours of equivalent widths for Hei and Hen lines using a model-grid with values for T e s and log g as indicated, and 
negligible wind densities, log Q — -14. Curves in grey color (left-hand side of Hei2.05 iso-contours) indicate net emission. Note 
that for the optical transition (Hei4471) the absorption increases as function of gra vity, whereas fo r the NIR lines this behaviour 
is mostly reversed. Asterisks denote the position of the calculated models (see also lPuls et alJ200^) . 
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Fig. 3. Line broadening profile functions (convolution of Doppler with Stark profiles, in [s -1 ]) for H r and BrlO, as a function 
of frequency displacement from the line centre in units of thermal Doppler-width. Both profiles have been calculated in Griem 
approximation, for an electron temperature T e = 40,000 K and three different electron densities typical for the line forming 
region. The grey line corresponds to a pure Doppler profile. Note that for H y the core of the Stark profile is identical with the 
pure Doppler profile, whereas for BrlO the Stark core is extremely density dependent and coincides with the Doppler profile only 
at lowest densities (see text). 
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lated for typical line-forming parameters, T e = 40,000 K and 
logn e = 11.5,12.5 and 13.5, respectively. The correspond- 
ing pure Doppler profile is overplotted (in grey). The decisive 
point is, that for H y , with relatively low upper principal quan- 
tum number, the Stark width is not considerably large, and the 
core of the profile is dominated by Doppler-broadening, inde- 
pendent of electron density. Only in the far wings does the well 
known dependence on n e become visible. On the other hand, 
for BrlO the Stark width becomes substantial (being propor- 
tional to the fourth power of upper principal quantum num- 
ber), and even the Stark-core becomes extremely density de- 
pendent. Only at lowest densities, the profile coincides with 
the pure Doppler profile, whereas for larger densities the pro- 
file function (and thus the frequential line opacity) decreases 
with increasing density. In the far cores, finally, the conven- 
tional result (<f>(v) correlated with n e ) is recovered. Thus, as a 
consequence of the dependence of Stark-broadening on den- 
sity, the line cores of the hydrogen lines with large upper prin- 
cipal quantum number become weaker with increasing gravity. 
Br y (with upper quantum number n=7) is less sensitive to this 
effect, cf. Fig^ 

In Fig. we demonstrate the different reactions of the Stark 
profiles on electron density (gravity) by comparing the syn- 
thesized emergent profiles of the high- and low-gravity model 
at r e ff=40,000 K, as described above. In particular, we com- 
pare these profiles with the corresponding profiles calculated 
with pure Doppler-broadening. For H r , the core of the Stark- 
broadened profile agrees well with the Doppler-broadened one 
(dashed) in both cases. The major difference is found in the 
far wings, which become wider and deeper as a function of 
(electron-) density, thus, being the most useful indicators for 
the effective gravity. For BrlO, on the other hand, the pure 
Doppler profile is much deeper than the Stark-broadened core, 
where the differences are more pronounced for the high gravity 
model. Note particularly that the (absolute) e.w. is larger for the 
low gravity model (although the high gravity model has more 
extended wings), since the major part of the profile is domi- 
nated by the core which is deeper for lower gravities! 

Actually, the same effect is already visible in the optical, 
namely for the prominent Hen lines at 4200 A (transition 4-11) 
and 4541A (transition 4-9, not shown here). The increase in 
absolute e.w. as a function of gravity is solely due to the wings. 
In accordance with BrlO/Brl 1, however, the cores of the lines 
become shallower with increasing gravity, not because of an 
effect of less absorbers, but because of less frequential opacity 
due to a strongly decreased broadening function. 

Let us allude to an interesting by-product of our investi- 
gation. A comparison of our synthetic NIR profiles with the 
observations will show that in a number of cases the observed 
BrlO/Brl 1 profiles cannot be fitted in parallel. In this case the 
line formation is well understood and the profiles from cmfgen 
are ident ical (note that also th e optical hydrogen lines agree 
well, see Rep olust et al]l2004T) . giving us confidence that our 
occupation numbers are reasonable and that the obvious differ- 
ences are due to inadequate broadening functions. 

On the other hand, since Hen4200/4541 is affected by al- 
most identical line broadening, we would like to suggest a so- 
lution for a long standing problem in the optical spectroscopy 
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Fig. 4. Influence of Stark-broadening for lines with low- and 
high-lying upper level as a function of gravity. Upper level, 
solid lines: Synthetic spectra of H r , BrlO and Hen4200 (n=4 
— >1 1) for an atmospheric model with r e ff= 40,000 K, log g=4.5 
and log Q = -14 (cf. dashed lines in Fig. [Q. Overplotted 
(dashed) are the corresponding profiles with pure Doppler- 
broadening. Lower panel: As the upper panel, but for an at- 
mosphere with logg=3.7 (solid lines in Fig.Q. For compari- 
son, the results for the Stark-broadened profiles from the upper 
panel are overplotted in grey. Note the similarity in effects be- 
tween BrlO and Hen4200. 



of hot stars: It is well known that for a wide range of O- 
star parameters the theoretical simulations of these lines (by 
means of both plane-parallel and extended atmospheres) have 
never been able to r eproduce the observations in parallel (e.g., 
Herre ro et al .120021) . where the largest discrepancies have been 
found in the line cores. The origin of this discrepancy is still 
unknown. 3 Due to the similarity of this problem to the one 
shown by Brl0/ll and accounting for the similar physics, we 
suggest that also in this case we suffer from an insufficient de- 
scription of presently available broadenin g functions (w hich 
are de scribed within the VCS-approach, see lSchoning & ButleJ 
1989j). Thus, a re-investigation of line broadening for transi- 
tions with high lying upper levels seems to be urgently re- 
quired. 

In summary, due to their tight coupling with electron den- 
sity, the cores of Brl0/ll and Henl.69/2.18 are excellent indi- 
cators of gravity, where deeper cores indicate lower gravities 
(if the (projected) rotational velocities are similar). 



3 Note that this problem is most probably not related to the presence 
of the Nm blend in Hen4200, since it occurs also in hot objects. 
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Fig. 5. Left panel: NLTE-departure coefficients for the lower (solid) and upper (dotted) levels of Heil.70 for a model with T e g= 
30,000 K, logg=3.4 and negligible wind. Overplotted in grey are the corresponding values for a similar model, but with lower 
gravity, logg=3.0. 

Right panel: As the left panel, but for the corresponding line source functions in units of the emergent continuum at 1.70/im. 



4.2. Hei lines: Influence of NLTE effects 

The peculiar behaviour of the hydrogenic lines could be traced 
down to the influence of the profile-functions, whereas the for- 
mation of most of the NIR Hei lines is dominated by NLTE - 
effects. As has been extensively discussed by Miha lasl ill 978b. 
iKudritzkil Jl979h. iNaiarro et alJ dl998h . IPrzvbilla & Butlerl 
J2004 and Lenorzer et alJJ2004l) . the low value of hv/kT leads 
to the fact that even small departures from LTE become sub- 
stantially amplified in the IR (in contrast to the situation in the 
UV and optical). A typical example is given by the behaviour 
of Heil.70 at temperatures below r e fr« 35,000 K, cf. Fig. [2] 
(note, that cmfgen gives identical predictions). Again, this line 
becomes stronger for lower gravity, in contrast to the well 
known behaviour of optical lines (compare with the Hei4471 
iso-contours). 

Fig. 13 gives a first explanation, by means of two atmo- 
spheric models with T e g =30,000 K, logg=3.4 and 3.0, respec- 
tively, and (almost) no wind. The departure coefficient of the 
upper level, b u , of this transition (4d 3 D) is independent of 
gravity and has, in the line forming region, a value of roughly 
unity (strong coupling to the Hen ground-state). The lower 
level of this transition (3p 3 P°) is quite sensitive to the dif- 
ferent densities, i.e., being stronger overpopulated in the low- 
gravity model. Consequently, the line source function, being 
roughly proportional to b u /b/, is considerably lower through- 
out the photosphere (right panel of Fig.|5|i, and thus the profile 
is deeper, even if the formation depth is reached at larger values 

of Tr oss . 

The reason for this stronger overpopulation at lower log g- 
values is explained by considering the most important pro- 
cesses which populate the 3p-level. First, the influence of col- 
lisions is larger at higher densities, which drives the departure 
coefficient into LTE. Second, the level is strongly coupled to 
the triplet "ground state" (i.e., the lowest meta-stable state) 
which, in the photosphere, is overpopulated as an inverse func- 
tion of the predominant density. The overpopulation (with the 



consequence of over-populating the 3/>level) is triggered by 
the strength of the corresponding ionizing fluxes. These are 
located in the near UV (roughly at 2,600 A) and are larger 
for high gravity models than for low gravity ones. This is be- 
cause the stronger Lyman-jump and the stronger EUV flux- 
blocking (higher densities — > lower metal ionization stages — > 
more lines) have to be compensated for on the red side of the 
flux-maximum to achieve flux conservation. 

If the ionization/recombination rates are dominating, the 
(photospheric) departure coefficients inversely scale with the 
flux at th e corresp onding edge (for similar electron tempera- 
tures, cf. lMihalaslll978l) . and for higher gravities we obtain 
lower departure coefficients (more ionization) than for lower 
gravities. Thus, the increase of the Heil.70 line flux with grav- 
ity is a final consequence of the different near UV radiation 
temperatures as a function of gravity. 

One might wonder why the strength of Hei447 1 is "well" 
behaved, since this line has the same upper level as Heil.70, 
and the lower level (2p 3 P°) is strongly coupled to the triplet 
ground-state as well. Actually, a simple simulation shows that 
for this transition the same effect as for the Hei 1.70 line would 
be present if the transition were situated in the IR. Only be- 
cause the transition is located in the optical (hv/kT » 1), the 
corresponding source functions are much less dependent on 
gravity (the non-linear response discussed above is largely sup- 
pressed). The profiles react almost only on the opacity, which 
is lower for lower gravity due to the lower number of available 
Hei ions. 

In summary, the Hei line formation in the optical is primar- 
ily controlled by different formation regions, since the source 
functions do not strongly depend on gravity, whereas in the IR 
the deviations from LTE become decisive. In particular, the in- 
fluence of considerably different source-functions is stronger 
than the different formation depths, where these source func- 
tions are larger for high-gravity models due to a less overpop- 
ulated lower level. 
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With respect to the singlet transitions (Hei2.05, react- 
ing inversely to the red compon ent of Hei2.l l ), we refer 
the reader to the d iscussion by iNaiarro et alJ i 19941) and 
Lenorzer et al . (2004). But we would like to mention that for 
a large range of parameters Hei2.05 reacts similar to the way 
described above, simply because the ionization/recombination 
rates (over-)populating the lower level (2s 'S, again a meta- 
stable level, located at roughly 3,100 A) remain the decisive 
ingredients controlling the corresponding source functions. 

Although the basic reaction of Hei2.05 on gravity is readily 
understood fr om this argum ent, we like to point out the follow- 
ing fcf. lNaiarroetalll9 94). For most of the objects considered 
here, the upper level of this transition (2p l P°) is populated 
by the Hei EUV resonance line at 584 A and other resonance 
lines (subsequently decaying to this level). Thus, it strongly 
depends on a correct description of line-blocking in this wave- 
length range, particularly on lines overlapping with the res- 
onance transitions. Moreover, any effect modifying the EUV 
will have a large impact on Hei2.05, e.g., wind clumping, if 
present close to the photosphere. This makes the line generally 
risky for classification purposes and the determination of T e g 
in those stars where Hen is no longer present (cf. Sect. 16. 1> . 

5. Comparison with results by Lenorzer et al. 

As already mentioned, iLenorzer et al.l (|Jo04) recently pre- 
sented a first calibration of the spectral properties of normal 
OB stars using near infrared lines. The analysis was based 
on a grid of 30 line-blanketed unified atmospheres computed 
with cmfgen. They presented 10 models per luminosity class 
I, III, a nd V, where wind -properties according to the predic- 
tions bv lVink et alJ |2000 ) have been used, and T e g ranges from 
24,000 K up to 49,000 K (cf. Fig.Q. Emphasis was put on the 
behaviour of the equivalent widths of the 20 strongest lines of 
H and He in the J, H, K and L band. For detailed information 
on procedure and results see lLenorzer et alJ ll2004l) . In order to 
check our results obtained by means of fastwind, we have cal- 
culated models with identical parameters and s ynthesized the 
same set of NIR lines (see also lPuls et all2 005). Note that cm- 
fgen uses a constant photospheric scale height (in contrast to 
fastwind), so that the photospheric density structures are some- 
what different, particularly for low gravity models where the 
influence of the photospheric line pressure becomes decisive. 

Since Lenorzer et al. calculated their hydrogenic profiles 
with th e erroneous broadening functions provided by Lemke 
ill997l) . the H and K band profiles have been recalculated by 
means of the Griem approximation by one of us (R.M.). The 
differences in the equivalent widths for the dwarf grid are 
shown in Fig. [6] In all cases the equivalent widths became 
larger, mainly because of increased line wings, and the most 
significant changes occurred for Brl 1 at lower temperatures. 
Note, however, that also Br r has become stronger throughout 
the complete grid. 

In Fig.0we now compare the detailed profiles of the strate- 
gic lines located in the H and K band of the present investi- 
gation (results from cmfgen in grey). The agreement between 
the results for the (almost purely photospheric) lines in the H 
band (BiTO/11, Heil.70 and Henl.69) is nearly perfect. The 
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Fig. 6. Comparison of equivalent widths (defined here in the 
conventional way) fo r NIR hydrogen lines of O-type dwarfs 
jLenorzer et alJ l2004|) using the (erroneous) broadening func- 
tions bv lL^rrAa ( [l997l) with results using broadening functions 
in the Griem approximation (dashed). Squares, circles and tri- 
angles correspond to Brll, BrlO and Br r , respectively. Note 
that in all cases the equivalent width becomes larger and that 
the most significant differences occur for Brl 1 at lower tem- 
peratures. 



only differences occur in the cores of BrlO, where cmfgen pre- 
dicts some emission for hotter objects, and some marginal dif- 
ferences in the far wings of the supergiants, which we attribute 
to the somewhat different density stratification in the photo- 
sphere. Additionally, cmfgen predicts slightly stronger Henl. 69 
lines for the hottest objects (models "IV" and "Ha") and for the 
supergiant model "61a". 

Concerning the K band, the comparison is also rather sat- 
isfactory, except for Hei2.05 at intermediate spectral type (see 
below). Concerning Br r , the dwarf models give rather similar 
results, with the exception of intermediate spectral types, where 
fastwind produces some central emission. We have convinced 
ourselves that this prediction is very stable (and not depending 
on any temperature inversion), resulting from some intermedi- 
ate layers where the population of the hydrogen levels is similar 
to the nebular case. Here the departure coefficients of the indi- 
vidual levels increase as a function of quantum number. In such 
a situation, a central emission owing to a strong source function 
is inevitable. For the supergiants, the major differences regard, 
again, the line cores, with cmfgen predicting more refilling. 

Somewhat larger differences are found for Hen2.18, again 
(cf. Henl.69) for the hottest models where cmfgen predicts sig- 
nificantly more absorption. 

Concerning Hei, the situation for the triplet line (blue com- 
ponent of Hei2.1 1) is as perfect as for Heil.70. The differences 
for the Hei (triplet) component located at the blue of Br y are 
quite interesting, cmfgen predicts an emission for hot stars but 
"nothing" for cooler objects, whereas fastwind predicts a rather 
strong absorption at cooler temperatures. To our present knowl- 
edge, this is the only discrepancy we have found so far (includ- 
ing the optical range) for a triplet line. 
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Fig. 7 . Comparison o f synthetic NIR lines for the grid of O-type dwarfs (upper panel) and supergiants (lower panel) described by 
Lenorzer et al .(2004), as calculated by fastwind and cmfgen (in grey). For hydrogen and Hen, the results reported by Lenorzer et 
al. have been recalculated using the Griem approximation. The horizontal and vertical lines in the bottom right corner indicate the 
scale used and correspond to 0.01/im in wavelength and 0.1 in units of the continuum, respectively. To simplify the comparison, 
the synthetic profiles have been convolved with a rotational profile corresponding to v sin i = 80 km s~'and degraded to a typical 
resolution of 10,000. 
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The only important discrepancies concern the Hei singlets 
(Hei2.05 and the red component of Hei2.11) of the super- 
giant and dwarf models in the range between models 5 to 
7. Starting from the hotter side, cmfgen predicts strong ab- 
sorption, which abruptly switches into emission at models 
no. 7, whereas fastwind predicts a smooth transition from 
strong absorption at model 4 to strong emission at model 8. 
Reassuring is the fact that at least the inverse behaviour be- 
tween H ei2.05 and Hei2.1 l(re d) (as discussed in lNaiarro et alJ 
1 19941 and lLenorzer et alJ .2004 ) is always present. 

Analogous c omparisons performed in the optical 
fchils et al] l2005t) have revealed that the strongest dis- 
crepancies are found in the same range of spectral types. The 
triplets agree perfectly, whereas the singlets disagree, because 
they are predicted to be much shallower by cmfgen than 
the ones resulting from fastwind. Again, the transition from 
shallow to deep profiles (at late spectral type) occurs abruptly 

in CMFGEN. 

IPuls et all (l2005t) have discussed a number of possibilities 
which might be responsible for the obvious discrepancy, but at 
present the situation remains unclear. One might speculate that 
this difference is due to subtle differences in the EUV, affecting 
the Hei resonance lines and thus a number of singlet states, as 
outlined in the previous section. We will, of course, continue in 
our effort to clarify this inconsistency. 

In addition to the detailed comparison performed in the H 
and band, we have also compared the resulting e.w.'s of some 
other important lines in the J and L band. Most important is the 
comparison for Br (y , w hich is a primary indi cator of mass-loss, 
as already discussed in lLenorzer et alJ l)2004h . Fig.[8]compares 
the corresponding e.w.'s, as a function of "equivalent width in- 
variant" Q' (see Lenorzer et al.). Generally, the comparison is 
satisfactory, and particularly the differences at large mass-loss 
rates are not worrying, since in this range the net-emission re- 
acts strongly on small changes in M. Real differences are found 
only for the weakest winds, probably related to "uncorrected" 
broadening functions used by Lenorzer et al. 

As we have found for Hen2.18, also the differences for the 
other Hen lines are significant. Note that we can only compare 
the e.w.'s and that the broadening as calculated by Lenorzer et 
al. suffers from erro neous line broadening . For Hen5-7 the be- 
haviour compared to lLenorzer et alJ 112004) is the same, but our 
lines are twice as strong in absorption in the case of giants and 
dwarfs. For the supergiants we obtain the strongest absorption 
at 36 kK, in contrast to 42 kK in the comparison models. The 
only difference found for Hen6-ll concerns the behaviour of 
supergiant and dwarf line trends. The supergiants in the com- 
parison models show stronger absorption lines than the dwarfs, 
whereas in our case the situation is reversed. 

At high r e ff, the models for Hen7-13 display a monotonic 
behaviour with the hottest models showing the strongest ab- 
sorption. Our hottest models display weaker absorption profiles 
(as was found in the detailed comparison of Hen 1.69 (7-12)), 
partly due to emission in the line cores. Finally, our emission 
lines obtained for Hen6-7 are twi ce as strong in the cas e of su- 
pergiants and giants compared to lLenorzer et alJ J2004I) . 

In summary we conclude that at least from a theoretical 
point of view, all H/K band lines synthesized by fastwind can 
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be trusted, except for Hei2.05 at intermediate spectral type 
and maybe Hen2.18, where certain discrepancies are found 
in comparison with cmfgen, mostly at hottest temperatures. 
Concerning the discrepancies of Hen in other bands, we have to 
clarify the influence of correct broadening functions, whereas 
for the Hei singlet problem work is already in progress. 

6. Analysis 

6.1. General remarks 

It might be questioned to what extent all decisive stellar and 
wind parameters can be obtained from a lone IR-analysis in 
the H and K band. In view of the available number of strate- 
gic lines, however, in most cases we are able to obtain the full 
parameter set, except for 

i) the terminal velocity, which in most cases cannot be de- 
rived from the optical either, and has been taken from UV- 
measurements. For our analysis, we have used the values 
given in the publications corresponding to sample I to III. 
The t erminal veloci t ies of sample IV have been adopted 
from Howart h et all Jl997h . If no information is (or will 
be) present, calibrati ons of Vn, as a function o f spectral type 
have to be used, e.g JKudritzki & Pulj(l2000l) . 

ii) the stellar radius, which can be infe rred from M v and the 
theoretical fluxes (Kudritzki 1980), and has been taken 
from the optical analyses in the present work. In future in- 
vestigations when no optical data will be available, a sim- 
ilar strategy exploiting infrared colors can certainly be es- 
tablished. 

In particular, BiTO and Brl 1 give clues on the gravity (if T e g is 
known), Hei and Hen define temperature and helium content, 
and Br y can serve as an M indicator, at least in principle. In 
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those cases, where only one ionization stage of helium is vis- 
ible, the determination of Fee becomes problematic, and also 
the uncertainty for T e g increases (see below). Due to the high 
quality of our spectra, however, both Hen lines are visible for 
most spectral types. 

Only for the coolest objects Hen vanishes, which occurs 
for spectral types later than 09 for dwarfs, about BO for giants, 
and again BO for supergiants (cf. Figs. l9lto ll li . In those cases 
it still should be possible to derive (somewhat more inaccu- 
rate) estimates for T e g, at least if some guess for Yu e is present. 
This possibility is due to the behaviour of Heil.70 vs. Hei2.05 
(Fig. |3, since the former line is almost only dependent on log g, 
whereas the latter depends strongly on r e ff (with all the caveats 
given earlier on). Unfortunately, the data for Hei2.05 are not of 
sufficient quality (except for HD 190864 and t Sco, where the 
latter just lies in the critical domain) that we could exploit this 
behaviour only once and had to refrain from an analysis of the 
remaining coolest objects (four in total). 

Because of the independence of Heil.70 on T e ff and the 
fact that BiTO/1 1 can always be fitted for certain combinations 
of r e ff/log g, a perfect fit in combination with completely erro- 
neous parameters would result if Hei2.05 had to be discarded. 
This is indicated in Fig.^2f° r HD 141 34, being a B3Ia super - 
giant (with r eff roughly at 18,000 K, see lKudritzki et all 19 99). 
which could be fitted with r e ff = 25,000 K. If, on the other 
hand, Hei2.05 had been available, the appropriate parameters 
should have been obtained, at least when the helium content 
could have been guessed. Such a guess of the helium abun- 
dance should always be possible for objects we are eventually 
aiming at in our project (cf. Sect. Q), i.e., for very young, un- 
evolved stars with unprocessed helium. 

Micro-turbulence. In agreement with the findings by 
IRepolust et alJ ll2004h . we have adopted a micro-turbulence 
of Vmrb = 10km s _1 for all stars with spectral type 07 or 
later regardless of their luminosity class, whereas for hotter 
O-type stars the micro-turbulent velocity has almost no effect 
on the analysis and we have neglected it. At spectral type 
06.5, our IR-analysis of HD 190864 (06.5 III) indicated that 
a micro-turbulence is still needed, whereas from 07 onwards 
Vturb did not play any role, e.g., for HD 192639 (07 lb). Since 
the former and the latter stars have r e ff = 37 and 35 kK, 
respectively, we conclude that at roughly r e (j = 36 kK the 
influence of v tur b on the H/He lines is vanishing, in agreement 
with our previous findings from the optical. 

Rotational velocities. For the (projected) rotational veloc- 
ities, _we_Jiave J _as a firs t guess, used t h e values pro vided 
bv IRepolust et all fc004 . iHerrero et al] feOOOl l2002fe and 
iHowarth et al.l i 19971) for sample I, II and III/IV, respectively. 
In our spirit to rely on IR data alone, we have subsequently in- 
ferred the rotational velocity from the (narrow) He lines, with 
most emphasis on Hei. Concerning sample I, the results from 
our IR-analysis are consistent with the velocities derived from 
the optical, except for HD 190864 and HD 192639, where the 
profiles indicated slightly lower values (10% and 20%, respec- 
tively), which have been used instead of the "original" ones. 



For sa mple II stars, in 3 out of 5 ca ses the "optical" values 
derived bv lHerreroetail d2000t [2002) were inconsistent with 
our IR-data. In particular, for HD 5689 we found a velocity of 
220 km s- 1 (instead of 250 km s" 1 ), for HD 15570 a veloc- 
ity of 120 km s" 1 (instead of 105 km s" 1 ) and for Cyg OB2#7 
our analysis produced the largest differences, namely V r sin i = 
145 km s~', compared to a value of 105 km s _1 provided by 
iHerreroet all d2002h (30% difference!). 

The values taken from IHowarth et al ] ill 9971) for the re- 
maining sample III/IV objects, finally, agree fairly well with 
our IR data, and are a lso consistent with the values derived by 
iKudritzki et alJ Jl999l) in their analysis of sample III objects. 

Let us finally mention that in those cases when Br y does not 
show emission wings, a statement concerning the velocity field 
parameter, is not possible, as is true for the optical anal- 
ysis. In order to allow for a meaningful comparison with re- 
spect to optical determinations of M, we have used the corre- 
sponding values derived or adopted from the optical. In future 
analyses, of course, this possibility will no longer be present, 
and we have to rely on our accumulated knowledge, i.e., we 
will have to adopt "reasonable" values for with all r elated 
problems concernin g the accuracy of M ("cf. IPuls et al.lll996t 
IMarkova et all2004 . 

6.2. Fitting strategy and line trends 

In order to obtain reliable fits, we applied the following strat- 
egy. At first, we searched for a coarse determination of the rele- 
vant sub-volume in parameter space by comparing the observed 
pro files with our large grid of synthetic profiles as described by 
IPuls et all d2005l) . which has a typical resolution of 2,500 K in 
T e fF, 0.3 in log g, and 0.25 in log Q. A subsequent fine fit is ob- 
tained by modifying the parameters by hand (using the "actual" 
values for R* and v M to obtain information on M additionally 
to log Q), where typically 10 trials are enough to provide a best 
compromise. In those cases, where at present no information 
about R„ is available (which concerns the three objects pre- 
sented in Table |3, "only" log Q can be derived. For the actual 
fits of these three objects we have, of course, used prototypical 
parameters for R* and v ro . Further discussion of related uncer- 
tainties is given in Sect.0 

Most weight has been given to the fits of the He lines 
(which are rather uncontaminated from errors in both broaden- 
ing functions and reduction of the observed material) followed 
by the photospheric hydrogen lines, BiTO/1 1, which sometimes 
strongly suffer from both defects. Least weight has been given 
to Br r , because of the nu mber of problem s inherent to this line, 
as recently des cribed by iLenorzer et alJ (|2004) and indepen- 
dently found bv l.Tokuthvlll2002l) . Particularly, the synthetic pro- 
files for larger wind densities, predicted by both fastwind and 
cmfgen, are of P Cygni type, whereas the observations show an 
almost pure emission profile. Moreover, from a comparison of 
equivalent widths, it has turned out that in a lot of cases the pre- 
dicted e.w. is much larger than the observed one, which would 
indicate that the models underestimate the wind-emission (re- 
member, that Br y forms inside the H a sphere). Often, however, 
this larger e.w. is due to the predicted P Cygni absorption com- 
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ponent which is missing in the observations, and we tried to 
concentrate on the Br r line wings in our fits ignoring any dis- 
crepancy concerning the predicted P Cygni troughs. If the syn- 
thetic lines actually predicted too few wind emission, this prob- 
lem would become severe for lines where pure absorption lines 
are observed, and should lead to an overestimate of M. We will 
come back to this point in the discussion of our analysis. 

Another important point to make concerns the Hei A2.W 
line (comprising the Hei triplet /12.1120 and the Hei singlet 
/12.1132). Close to its central frequency, a broad emission fea- 
ture can be seen (at ^12. 115) in the spectra of hot stars. This 
line can either be identified as Nm (n = 7 — >8 ) or as Cm (n 
= 7 — >8) or maybe both ( lHanson et alJ 1 19961: iNaiarro et alJ 
11997a! 2004). 4 This feature is seen in stars of all luminosity 
classes, for stars hotter than and including spectral type 08 in 
the case of dwarfs and giants and 09 in the case of supergiants 
(though its designation is somewhat unclear, as Hei2. 1 1 resem- 
bles a P Cygni profile in late-O supergiants, possibly mimick- 
ing this feature). Since our present version of fastwind syn- 
thesizes "only" H/He lines and their analysis is the scope of 
the present paper, we are not able to fit this feature, but have 
to consider the fact that this feature significantly contaminates 
Hei2.11. 

Due to the well-resolved spectra, the two Hei lines over- 
lapping with Br y as mentioned in Sect. |2 i.e., the Hei triplet 
/12.1607 and the Hei singlet A2. 16 17, are also visible in certain 
domains. For supergiants later than 05, Hei2.1607 begins to 
appear in the blue wings of Br y , and in two stars, HD 30614 
and HD 37128, the Hei2.1617 singlet seems to be present, even 
if difficult to see. In the giant spectra, Hei2.1607 can be seen 
from spectral type 09 onwards, and in the dwarf spectra this 
line appears in spectral types later than 08. 

The strength of the Brackett lines in supergiants (Fig. II 11 
shows a smooth behaviour as a function of spectral type, apart 
from certain fluctuations such as blends in the late O- and early 
B-type stars. As one moves from early B-type to mid O-type 
(i.e., 05), the Br r absorption weakens, and from mid to the 
earliest O-types the line profiles switch into emission, where 
the emission at the blue wings of Br r is much more pronounced 
(except for HD 15570), presumably due to the overlapping Hen 
blend. 

As for the photospheric Brl0/1 1 lines, we can see that these 
absorption profiles show an extremely continuous behaviour, 
being rather weak for early O-type stars and increasing in 
strength towards early B-types. Hence, the cooler supergiants 
show the most prominent and sharpest absorption features. The 
emission features visible at the blue side of BrlO in the hottest 
supergiants are due to an unidentified feature. 

Fig. II llshows that the observed Brl0/1 1 profiles are mostly 
well reproduced by the theoretical predictions, although at hot- 

4 Due to the rather similar structure and the fact that these tran- 
sitions occur between high lying levels, the predicted transition fre- 
quencies are almost equal. Since most of the stars in the Olf phase 
will have depleted C and enhanced N, however, the major contribution 
should be due to Nm and possibly also due to Nv2.105 (10 — >11) for 
the hottest objects (F. Najarro, priv. comm.). Cm will be contributing 
if Civ at 2.07-2.08yum is strong. 



ter temperatures certain inconsistencies arise, particularly with 
respect to the line cores. Most interestingly, in a number of 
cases we could not fit both profiles in parallel, and typically 
Brl 1 is then of better quality. Since we have convinced our- 
selves that the differences most probably are not exclusively 
due to reduction problems, we repeat our hypothesis that the 
broadening functions are somewhat erroneous, cf. Sect. 14.11 
Again, for the theoretical profiles for Br r , we would like to 
mention that for emission lines the wings are fairly well repro- 
duced in contrast to the line cores. 

The Heil.70 line shows a very smooth behaviour, being ab- 
sent in the hottest and most luminous star, CygOB2 #7, and 
successively increasing towards late O-type and early B-type 
stars. This also applies to the sharpness of the profiles. As has 
been stressed earlier on, Henl.69 and Hen2.18 vanish in su- 
pergiants of spectral type B0 (being still detectable for a Cam, 
09.51a) 

The situation is similar in the case of giants (Fig. II Oi and 
dwarfs (Fig.|9). All hydrogen and Hei lines show the systematic 
variations expected, namely, an increase in strength from early 
O-types to early B-types. The model predictions do agree well 
with the observed profiles, again, except for certain discrepan- 
cies between BrlO vs. Br 11. Since Br r remains in absorption 
throughout the entire spectral range, it can be reasonably fitted 
in most cases (whether at the "correct" value, will be clarified 
in Sect. Particularly, the Hen profiles give almost perfect 
fits except for very few outliers, and vanish at 09 for dwarfs 
and about B0 for giants. 

6.3. Comments on the individual objects 

In the following, we will comment on the fits for the individ- 
ual objects where necessary. Further information on the objects 
can be found in the corresponding publications concerning the 
optical analyses, see Tab. ^ A summary of all derived values 
can be found in Tables|2]and|4] 

6.3.1. Dwarfs 

HD 64568. The fit quality of the lines is generally good ex- 
cept for Br r . The theoretical profile displays a central emission 
which is more due to an overpopulated upper level than due 
to wind effects (cf. Sect. [5} and, thus, cannot be removed by 
adopting a lower M. Moreover, the theoretical Hen line would 
become too strong if a lower M were used. Insofar, the present 
fits display the best compromise. Since no radius information 
is available, only the optical depth invariant log Q is presented. 

HD 46223 and HD46150. For both objects, the fit quality is 
satisfactory, except for the Hen lines (particularly in HD 46 1 50) 
and Br r . The former lines are predicted to be too strong and for 
the latter there is, again, too much central emission present. 
It is, however, not possible to reduce the temperature in order 
to fit the Hen line, since this would adversely affect Hei. A 
further reduction of Tee (adopted here to be "solar", i.e., 0.1) 
is implausible, so that the presented line fits reveal the best fit 
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Fig. 9. Line fits for hot dwarfs with spectral types ranging from 03 to 07 (upper panel) and cool dwarfs with spectral types 
ranging from 08 to B2 (lower panel). The lowermost object (HD 36166, B2V) has not been analyzed due to missing Hen and 
Hei2.05 (see text). The horizontal and vertical lines in the bottom right corner indicate the scale used and correspond to 0.01 
microns in wavelength and 0. 10 in units of the continuum, respectively. 
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quality possible. For HD 46223 we can derive only the optical 
depth invariant log Q due to missing radius information. 



HD 1 5629. The fit quality for the He lines is very good, and we 
confirm the helium deficiency to be Fh r = 0.08 as determined in 
the optical, see Repolus t et alJ (12004). Br r , again, suffers from 
too much central emission, and the cores of Brl0/1 1 are much 
narrower than predicted (at least partly, as some of the narrow- 
ness might be due to reduction problems). The mass-loss rate 
is moderate with a value of M » 1.3 T0~ 6 M o /yr which repre- 
sents the same value as determined in the optical, whereas log g 
is found to be larger by 0. 1 dex. 



HD5689. Again, moderate mismatches for the H lines are 
found, whereas the He lines provide a good fit. Br r does not 
show a central emission anymore, but the theoretical profile 
seems to be too broad. The same problem (very steep increase 
on the blue side, almost perfect fit on the red side) seems 
to be present also in HD 217086 (and, to a lesser extent, in 
HD 203064 and HD 191423), and we attribute some of this dis- 
agreement to reduction errors, although an underestimate of the 
Hen blend (which is in emission in this parameter range) might 



be possible as well. Since all four stars are very fast rotators, 
effects from differe ntial rotation in combinatio n with a non- 
spherical wind (cf. IPuls et all 1 19961 iReoolust et alJEoO^I and 
references therein) cannot be excluded, see below. 

In the case of Br 10/1 1, on the other hand, problems in the 
broadening functions might explain the disagreement, as al- 
ready discussed. Finally, the absorption trough of the theoret- 
ical profile for Hei2.11 seems to be too strong, but might be 
contaminated by the bluewards Nm/Ciii complex. 

HD 21 7086. A very similar fit quality as found for HD5689 
has been obtained for this star, although BiTO/11 are now in 
better agreement. The parameters determined are comparable 
to the ones obtained from the optical, including the overabun- 
dance of He (Fee = 0. 15). An upper limit for the mass-loss rate 
has been derived, which is less than half the value obtained 
from the optical. 

H D 1 3268. The theoretical prediction reproduces the observa- 
tion quite well, especially in the case of Heil.70 and both Hen 
lines. As for the hydrogen lines, the two photospheric lines 
BrlO and Brll show too much absorption in the line cores, 
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whereas Hei2.1 1 shows the same trend as already discussed for 
HD 5689. The fit quality for Br y , however, is much better, and 
even the Hei2. 1607 (triplet) blend is reasonably reproduced, al- 
though slightly too strong. For the mass-loss rate only an upper 
limit of 0.17 TO~ 6 M /yr can be given, for an adopted value of 
yS = 0.80. The enhanced helium abundance Yu e - 0.25, as found 
in the optical, could be confirmed, giving the best compromise 
regarding all He lines. 

HD 149757 and HD 37468. The very good fit quality makes 
further comments unnecessary. 

HD 1 49438. r Sco is probably one of the most interesting stars 
of the sample analyzed, since it is a very slow rotator and all 
features become visible at the obtained resolution. Although 
only the K band observation is available, it can be seen that we 
obtain a very good fit quality for all H and He lines present 
(Hen is absent at these temperatures). As discussed before, in 
those cases where only one ionization stage of helium is vis- 
ible, the determination of Fee becomes problematic, and also 
the uncertainty for T e ff increases. Since in the case of r Sco 
we could make use of the Hei2.05 line, we could still deter- 
mine the effective temperature (resulting in a similar value as 
in the optical), on the basis of an adopted value Yu e =0.1. Also 
the mass-loss rate is well constrained from the resolved central 
emission feature in Br r , having a value of 0.02 - l O^M o/vr. 
From a similar investigation by Przvb illa & Butlerl (|2004), ex- 
ploiting the central emissions of Pf r , Pfg and Br a as well, they 
derived a value of 0.009 •lO _6 M /yr (factor two lower) as a 
compromise, but have adopted a different velocity-field expo- 
nent (J3 = 2.4 instead of/? = 1.0 used here) and utilized the 
"canonical" value for \ogg = 4.25 which fits H y . In our case, 
however, and in the spirit to rely on a lone IR analysis, we pre- 
ferred a lower value, \ogg = 4.0, since in this case the emission 
feature is better reproduced (much narrower) than for a higher 
gravity, whereas the differences in Hei2.05 (and concerning the 
line wings of Br y !) are almost negligible. If we have had the in- 
formation on Brl0/1 1, this dichotomy could have been solved. 

Having finished our investigation, one of us (R.M.) has 
analyzed the optical spectrum of r Sco, also by means of 
fast w ind. Details will be published elsewhere (Mo kiem et alJ 
1200.4 . Most interestingly, he obtained perfect line fits, at pa- 
rameters T eS = 31,900 K, logg = 4.15, F He = 0.12 and M = 
0.02 . . . 0.06 -10~ 6 M o /yr (for velocity exponents = 2.4 ... 
0.8, respectively). We like to stress that this analysis has not 
been biased by our present results from the IR, since it was 
performed by an "automatic" line fitting procedure based on a 
genetic algorithm. 

HD 361 66. This object has not been analyzed, due to missing 
Hen and Hei2. 05 lines. 

6.3.2. Giants 

HD 1 5558. Also for this star, only the K band observation is 
available, and because of the high temperature and rather large 



M no independent information concerning r e jf and logg can 
be derived. Thus, we adopted the effective temperature at its 
"optical" value, T e ff = 41,000 K. With this value, a simultane- 
ous "fit" of log g, Fee and M resulted in the synthetic spectrum 
displayed. M was constrained by the wings of Br y , and Fn e = 
0.08 derived on the basis that at this value Hen is still some- 
what too strong, log g is rather badly defined, since a variation 
by +0.2 dex gives only small differences in all three observed 
lines. In conclusion, the fit obtained allows to reliably constrain 
the mass-loss rate alone, and this only the temperature actu- 
ally has the adopted value. Note, however, that a (much) lower 
value is excluded since the predicted Hen2.18 line would be- 
come too weak (cf. Fig|2] lower right panel). 

HD 1 90864. The analysis gives a consistent fit for all lines (in- 
cluding Hei2.05!) except for Br y , where the theoretical profile 
of Br y shows too much central emission. The parameters re- 
mained almost the same compared to the optical except for the 
helium abundance, Fe e , which has been increased from 0.15 to 
0.20. 

HD203064 and HD 191423. The analysis for HD 203064 
yields a consistent fit for all lines, except for Hei2.1 1 which dis- 
plays a similar problem as described for HD5689. We recov- 
ered the same values for r e jf and log g as in the optical, though 
the helium abundance had to be doubled and also the mass-loss 
rate increased by roughly 80%. The theoretical profile of Br y 
for both stars is slightly broader than observed, although the ef- 
fect is weaker than found for HD 5689 and HD 217086. Note in 
particular that for both giants H a turned out to be narrower than 
predicted, with "e mission humps" pres ent on both sides of the 
absorption trough ( Rep olust et al.lEo04l Fig. 6). Summarizing 
and considering their extreme rotation velocities (V r sin i being 
300 km s 'and 400 km s^for HD 203064 and HD 191423, re- 
spectively), our above hypothesis of rotational distortion is the 
most probable solution for the apparent dilemma in these cases. 

Also for all other lines, HD 191423 behaves very similarly 
to HD 203064, although a better fit quality for Hei2. 1 1 is found, 
while Heil.70 has become worse (we aimed for a compromise 
between both lines). 

6.3.3. Supergiants 

Cyg OB2 #7. This star, being the hottest one in the sample, 
shows an enormous discrepancy in the Br y line, due to the ob- 
served central emission, which is not predicted by our simu- 
lations. It is the only star in our sample where we find the 
same problem in Hen2.18, i.e., where the theoretical predic- 
tions with respect to its morphology could not be confirmed. 
In order to determine a fairly reliable mass-loss rate, we have 
concentrated exclusively on the wings of Br y . The parameters 
derived agree with their values from the optical, except for the 
helium abundance. The determination of this quantity is prob- 
lematic d ue to missing Hei. I n contrast to the optical value, 
F He = 0.3 jHerrero et all2002l) . our best fit favoured F He =0.1, 
whereas simulations using the optical value have adversely af- 
fected the H lines. Moreover, to preserve the good fit quality of 
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Fig. 11. As Fig.|5J but for hot supergiants of spectral type 03 to 07 (upper panel) and cool supergiants (09 to B3, lower panel). 
The three lowermost objects (HD 13854, Bllab, HD 13866, B2Ib and HD 14134, B3Ia) have not been analyzed. The synthetic 
profiles overplotted for HD 14134 show a perfect fit for completely wrong parameters (7 , e ff= 25,000 K, logg= 2.7), indicating 
that a spectroscopic H/K band analysis is impossible if Hen and/or Hei2.05 are missing (see text). 
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Henl.69, we would have to lower T e g significantly if Fee = 0.3 
were the correct value. (Actually, a temperature already lower 
by 1,500 K compared to the optical has been used to achieve 
the displayed fit). Interestingly, a re-analysis of Cyg OB2 #7 in 
the optical performed by one of us (R.M.) resulted in a value 
just in between, namely Fee = 0.21 (at T e ff = 46,000 K). The 
emission on the blue side of BrlO is due to an unknown feature, 
as discussed in Sect. 16.21 

HD 6681 1 . The fit quality is generally good, except for Br y , 
which again shows much more central emission than predicted. 
The wings, on the other hand, could be well fitted and gave a 
mass-loss rate of 8.8 TO~ 6 M /yr, in agreement with the optical 
value. BrlO is contaminated by an unknown feature on the blue 
side, but to a lesser extent than in Cyg OB2 #7. 

HD 15570 and HD 14947 show very similar profiles, and 
could be reasonably well fitted. Note the prominent emission 
in Br r . This could not be matched, so we had to concentrate on 
the wings. In both cases Hen2.18 gives an additional constraint 
on M, since at higher values the (theoretical) wings would show 
too much emission. 

Cyg OB2 #8C and Cyg OB2 #8A. These stars, being of rather 
similar type and displaying rather similar profiles (with the 
noticeable difference of Heil.70, immediately indicating that 
8A is somewhat cooler than 8C), have been carefully ana- 
lyzed in the optical (and, again, re-analyzed by R.M.). From 
the optical, both stars have significantly different gravities (well 
constrained from the Balmer line wings), where object 8C 
with loss = 3.8 ha s a rather large gravity for its type, cf. 
iHerrero et al] d2002h . The values derived from the IR, on the 
other hand, are much closer to each other, namely 3.62 and 
3.41, respectively. 5 According to the observed shape of the 
profiles and their corresponding theoretical fits, a higher log g 
would lead to severe inconsistencies. Apart from gravity, how- 
ever, the other parameters derived are comparable to their op- 
tical counterparts, including the differences in M, although the 
fit quality of Br y is dissatisfying. 

HD 1 92639. For this star, we found a reasonable compromise 
concerning the fit quality of the lines present. We derived a 
logg value of 3.3 compared to 3.45 in the optical, because of 
the wings Brl0/1 1 (note the different degree of inconsistency in 
the lines cores!) and due to the shape of Hen2.18. With a value 
of logg = 3.45 Hen2.18 becomes even narrower, with a more 
pronounced P-Cygni type profile. The helium abundance was 
raised to 0.3 (from 0.2 in the optical) in order to fit the Hei and 
Hen lines appropriately in combination with the derived T e s. 
Also in this case, the observed Br y line shows a central emis- 
sion which could by no means be reproduced. The Hei2.1607 

5 Recently. Ide Becker et alj i2004) have identified object 8A as an 
06 1/05.5 III binary system, therefore the derived parameters remain 
doubtful. In our spirit to compare with optical analyses, we trea ted the 
system as a single star, in accordance with Herrero et al. (2002). 



triplet blend showing up in the theoretical prediction is not yet 
present in the observation. 

HD 21 0809. Part of the observed discrepancy in Br y might 
be attributed to int rinsic variations in t he notoriously variable 
wind of this star (Mark ovaet ail 2005). though it is also pos- 
sible that some (though not all) of the mismatch arises from 
errors in the removal of the Br y feature in the telluric stan- 
dard. Fortunately, the line wings could be fitted fairly well, 
resulting in a mass-loss rate of 5.80 TO~ 6 M /yr compared to 
5.30 •lO _6 M /yr in the optical. The major difficulty encoun- 
tered was to fit the Hei and Hen lines simultaneously. In fact, 
a decrease in T e g leads to an even more pronounced P-Cygni 
type profile for Hen2.18 for the given mass-loss rate, as was 
already true for HD 192639. We regard our solution as the best 
compromise possible, accounting for the fact that by a reduc- 
tion in T e ff we would also increase the apparent dilemma in 
Brl0/1 1 and the Hei component in Br y . The helium abundance 
was raised by 0.06 to 0.2 in order to find a compromise for the 
He lines. 

HD30614. For this star a very good fit quality was obtained 
making further comments unnecessary. 

HD 209975. The stellar profiles are fairly well reproduced and 
represent the best compromise possible. All hydrogen features 
predicted are a little too strong, with some contamination on the 
blue side of the profiles. The parameters obtained are compara- 
ble to the optical ones, except for logg, where we determined 
a smaller value (0. 15 dex). 

HD37128 (e Ori). Almost perfect fit. Let us only point out 
that the derived value for T e s represents an upper limit, since 
from this star onwards Hen is no longer present and Hei be- 
comes rather insensitive to T e g, so that without Hei2.05 further 
conclusions are almost impossible. 

HD 13854 and HD 13866 have not been analyzed, due to 
missing Hen and Hei2.05. 

HD 14134. As above. The "theoretical" spectrum displayed 
in Fig II II shows the insensitivity of the Heil.70 and Hei2.11 
lines to T e ff for temperatures below 30,000 K. Although a vir- 
tually perfect fit has been obtained, the synthetic model (T e s = 
25,000 K, logg = 2.70) is located far away from realisti c val- 
ues (r oughly at r eff = 1 8,000 K, log g = 2.20, cf. lKudritzki et alJ 



7. The "log 2-approach" 

Before we compare our results from the IR with optical data, let 
us briefly consider those objects where no optical information 
is available. Tab. [2] summarizes the corresponding parameters 
which constitute a "by-product" of our investigations. Because 
of the missing radius information, we quote the correspond- 
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Table 2. Adopted and derived stellar and wind parameters obtained from spectra in the infrared for sample IV objects not analyzed 
in the optical. Since no radius information is available, only the optical depth invariant, log Q, can be derived (Eq.^. T e ff in kK, 
V r sin i in km s _1 . Centrifugal correction (— > log g trlK ) assuming a typical radius. 



Star 


Sp.Type 


Teff 




log gtrue 




V r sin ( 


log<2 


P 


HD 64568 


03 V((f)) 


45.0 


3.85 


3.86 


0.10 


150 


- 13.00 


0.90 


HD 46223 


04 V((f)) 


42.0 


3.70 


3.71 


0.10 


100 


-12.70 


0.90 


HD 37468 


09.5 V 


30.0 


4.00 


4.00 


0.10 


80 


-14.10 


1.00 



Table 3. Different models for f Pup at identical log Q - 
-12.02. The first entry corresponds to our reference model 
(for additional parameters, cf. Table[4}. For models A to D we 
have varied R t , v M and M in such a way as to preserve the 
optical-depth invariant Q whereas the other parameters remain 
unchanged. R* in R , in km s'and M in lO~ 6 M /yr (see 
text and Fig. 1121 . 



Model 


R, 


Voo 


M 


log 2 


HD 66811 


19.4 


2250 


8.77 


-12.02 


HD 668 11- A 


15.0 


1900 


4.56 


-12.02 


HD66811-B 


15.0 


2550 


7.19 


-12.02 


HD66811-C 


25.0 


1900 


9.94 


-12.02 


HD66811-D 


25.0 


2550 


15.46 


-12.02 



ing values for the optical depth invariant, log Q, instead of the 
mass-loss rate M. 

Though we will no longer comment on these stars in the 
following, we would like to point out that all derived param- 
eters appear to be fairly reasonable, except for the gravity of 
HD 46223, which is rather low for a dwarf of spectral type 04. 

It might be questioned how reliable these parameters are, 
given the fact that, more precisely, log Q is actually a scaling- 
quantity for recombinatio n lines formed in the wind. As out- 
lined bv lPuls et alJ J2005h . however, this quantity is a suitable 
compromise concerning the scaling properties of other impor- 
tant physical variables, namely density (oc M/(/?^.Voo)) and the 
optical depth of resonance lines from major ions (scaling via 
M /(R-kV^J). Thus, it might be used as a general scaling invari- 
ant, though only in an average sense. Because of these different 
scaling properties, it is to be expected that models with identi- 
cal log Q but different combinations of R„, v x and M result in 
somewhat different profiles, and in the following we will ex- 
plore the corresponding uncertainties. By this investigation, we 
will also clarify in how far uncertainties in stellar radii (i.e., dis- 
tances) and terminal velocities (which will be present in future 
applications when analyzing the IR alone, cf. Sect. 16. 1> might 
influence the derived parameters, particularly log Q. 

A firs t impr ession of this problem has been given by 
IPuls et all J2QQ51 Fig. 18) who discussed the variations of the 
synthetic optical spectrum of a Cam when the stellar radius is 
modified, while keeping the terminal velocity and Q. The dif- 
ferences turned out to be marginal, except for the innermost 
cores of a variety of lines. H a (being the prototypical recombi- 
nation wind line), on the other hand, showed perfect agreement 
with the original profile. 

We have performed a corresponding analysis in the IR, ad- 
ditionally allowing for variations in Voo, by means of our best 
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Fig. 12. Variation of synthetic IR spectra for different combina- 
tions of R*, Voo and M but identical optical depth invariant Q. 
Bold: profiles for our reference model of £ Pup. Grey: maxi- 
mum und minumum profiles resulting from the models A to D 
of Tableland the reference model, see text. The scale in the 
lower left corner corresponds to 0.01 microns in wavelength 
and 1 % of the continuum, respectively. Note that in almost all 
cases the maximum/minimum variations are below this value. 
All original profiles have been convolved with a typical rota- 
tional speed of v sin i = 100 km s _1 . 

fitting model for f Pup, cf. Table 0] Our choice was motivated 
by the fact that this object is one of the few which shows a cer- 
tain degree of wind-emission in the wings of Br y and Hen2.18, 
thus displaying profiles which have a significant contribution 
from the wind. We have calculated four additional models, de- 
noted by A to D, where the modified stellar/wind-parameters 
(constrained by the requirement of preserving Q) are displayed 
in Table [3] In particular, we have varied R, = 19.4 R Q by 
roughly + 25% and Voo = 2250 km s _1 by + 300 km s _1 , con- 
sistent with the typical uncertainties in R* and Voo if taken from 
calibrations. Within the four models, the mass-loss rates (with 
a reference value of 8.5 T0~ 6 M o /yr) range from M = 4.6 to 
15.5 -lO- 6 M /yr. 

In Fig. ^] we compare the resulting variations of the cor- 
responding profiles with our reference solution from Fig. ^2 
Since the differences turned out to be rather small, we display, 
at each frequency point, the minimum and maximum normal- 
ized flux value with respect to all five spectra, in order to il- 
lustrate the maximum uncertainty due to the various parameter 
combinations. As it is evident, in almost all cases this maxi- 
mum uncertainty lies well below ± 1 % of the continuum flux 
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level, at least if v sin i is not very low. Accounting for the ad- 
ditional noise within the observation, we would have derived 
almost identical stellar parameters including log Q (with dif- 
ferences well below the typical errors as discussed in the next 
section) if a fit by any of the additional models had been per- 
formed. Thus, our hypothesis that log Q can be used as a global 
scaling invariant is justified indeed. 

By comparing the individual spectra from the four models 
with our reference spectrum, it turned out that the major source 
of disagreement results from models A and D, whereas mod- 
els B and C show almost identical profiles. The reason for this 
discrepancy/agreement is readily explained by noting that the 
run of p as a function of tr oss (which most importantly con- 
trols the emergent spectrum) is very similar for the reference 
model and models B and C, whereas for model A higher den- 
sities and for model D lower densities (as a function of tr oss ) 
are found, particularly in the transition region between pho- 
tosphere and wind. Consequently, the recombination rates in 
model A are somewhat larger, leading to deeper Hei profiles, 
and vice versa for model D. The resulting uncertainty with re- 
spect to r e ff , however, lies well below + 500 K. 

Note that all these investigations have been performed by 
keeping /3 at its reference value. If (3 cannot be constrained 
from the emission wings of wind lines, the uncertainty in log Q 
can b ecome severe, extending to a factor of two ("cf. lPuls et alJ 
1996). 



8. Comparison with optical data 

In the following, we now return to the most important ques- 
tion, to what extent a lone near IR analysis is able to recover 
the parameters from an analogous optical analysis. The corre- 
sponding data can be found in Tab. |4] 

Figs.^]and^^compare the results for T e ff, \ogg and M for 
all stars which have been analyzed in the optical. In each fig- 
ure, we have indicated error bars usually quoted for the corre- 
sponding "optical" quantity. In particular, the maximum errors 
for effective temperature are on the order of +5% (correspond- 
ing to +2, 000 K at r eff = 40,000 K 6 ) and typical errors for 
logg are +0.1. The "error bars" for the mass-loss rates, indi- 
cated as +0.2 dex, correspond to mean values attributed to M 
measurements from H ff . Note, however, that the actual preci- 
sion is an increasing function of M, being hi gher than 0.2 de x 
for low M and lower for larger values (e.g.. |Puls et al.lll996i) . 
Remember also that all our simulations (both in the IR and the 
optical) have been performed with un-clumped models, i.e., the 
derived mass-loss rates represent upper limits and may need to 
be corrected. 

From the three figures, it can immediately be seen that the 
majority of IR- values are in reasonable agreement with the cor- 
responding optical data. Most importantly, no obvious trend is 
visible, neither as a function of the parameter itself nor as a 
function of luminosity class ("lc"). A weak trend in T eB cannot 
be excluded though: From T eB = 35,000 K on, the IR data are 



6 This uncertainty is also consistent with the uncertainty related 
to the Hei singlet p roblem possibly affecting our optical analyses, 



cf. lPuls etal J <2005t Sect. 10). 



distributed more towards lower values (than derived from the 
optical). 

In the following, we will briefly discuss the outliers, i.e., 
those objects which are located beyond the indicated error bars 
and thus must be interpreted as severe mismatches. 

With respect to T e ff, three objects seem to be discrepant, 
at least at first glance. Concerning the moderate deviation of 
HD 37128 (e Ori, with an optical T eB = 27,500 K), remember 
that our IR-value is an upper limit only. In particular, we have 
compared our results wi th those derived from a recent, line- 
blanketed re-analysis by Urbaneia (2004), exploiting the ion- 
ization balance of silicon (lower entries in Tab. |4). The other 
two outliers, HD 46150 (lc V, appearing to be 3,000 K cooler) 
and HD 15570 (lc I, 4,000 K cooler), are unproblematic as 
well, since the corresponding optical analyses have been per- 
formed by unblanketed models, and the obtained differences 
are just of the expected order of blanketing ef fects at — 
42,000 K, i.e., roughly A7 W * -3. 500 K fe.g.. lMarfins et ail 
l2002HRepolust et al.ll20 04). New positions corrected for blan- 
keting effects of this amount have been indicated by arrows. 

Blanketing effects do not only affect effective temperatures, 
but also gravities. In particular, one must still be able to fit 
the Balmer line wings, which, for the correspondingly cooler 
temperatures, usually yields lower values for \ogg. Because of 
possible differences in the derived He-abundance, however, the 
actual correction depends strongly on the specific situation, and 
corrections towards significantly lower v alues and of n egligible 
amount have been found in parallel, see Repo lust et alT(l2004 
Fig. 16). 

In order to obtain at least an impression about the situation, 
we have used the log g - T eB calibrations provided by Repolust 
et al. (their Fig. 17, see also Mar kova et al.ll2004 L These cali- 
brations have been derived from their analysis of Galactic O- 
stars, utilizing the same code as applied here (augmented by 
the results for O-dwarfs from M artins et al.ll2002l) . and should 
be valid within ±0.1 dex, at least differentially. Using the cor- 
rected values for T e g from above (i.e, 3,500 K cooler than the 
unblanketed ones), we find log g = 3.62 and 3.87 for HD 15570 
and HD46150, respectively, where the former value could be 
even lower, due to the extreme character of the object (04If+). 
The corrected positions have been indicated in Fig. [O] right 
panel, again by arrows. Obviously, the IR-gravity of HD 15550 
is now almost consistent with the calibrated value from the op- 
tical. For HD 46150, on the other hand, the situation has be- 
come worse, with a discrepancy of -0.15 dex, if our calibration 
applies for this star. Note, however, that also the unblanketed 
optical gravity is (very) low for a dwarf of this spectral type, 
which might indicate that this object is somewhat peculiar, and 
an optical re-analysis is certainly required to obtain firm con- 
clusions. 

t Sco has not been included into this comparison, since no 
photospheric hydrogen lines have been observed and its gravity 
is almost unconstrained with respect to our NIR analysis (cf. 
Sect. l6.3.TY Remember that we have favoured the lower value 
alone because of the shape of the central emission in Br y , which 
constitutes the only difference between a logg = 4.0 and \ogg 
= 4.25 model, given the observed lines. In future analyses with 
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Table 4. Comparison of stellar and wind parameters in the optical and the near infrared derived using fastwind. T e g in kK, /?» 
in R , M in 10~ 6 M o /y r. loss values a r e corre cted for centrifugal acceleration. If not explicitely indicated, the optical parameters 
have been taken from Repolust et al . (2004). Note that Cyg OB2 #8 A has been recently identified as an Q6 I/Q5.5 III binary 
system. In order to compare with previous analyses, we have retained its single star designation as used bv lHerrero etal.N2002l) 



Qtnr 


S p. Type 


D 
/V# 


7 eff 


T"ir 




l0 g £lrue 


v opt 
^He 


^He 


IVl r 


IVl 


Cyg OB2#7" 


03 If 


14.6 


45.5 


44.0 


3.71 


3.71 


0.2 1 D) -0.30 


0.10 


9.86 


10.00 


HD 66811 


04 I(n)f 


19.4 


39.0 


39.0 


3.59 


3.59 


0.20 


0.17 


8.80 


8.77 


HD 15570 2 ' 


04If+ 


22.0 


42.0 


38.0 


3.81 


3.51 


0.18 


0.15 


17.8 


15.20 


Cyg OB2#8C n 


05If 


13.3 


41.0 


39.0 


3.81 


3.62 


0.09 


0.10 


2.25 


2.00 


HD 14947 


05 If+ 


16.8 


37.5 


37.5 


3.48 


3.48 


0.20 


0.20 


8.52 


7.46 


Cyg OB2 #SA [} 


05.5 1(f) 


27.9 


38.5 


37.0 


3.51 


3.41 


0.10 


0.10 


13.5 


11.50 


HD 192639 


07 lb 


18.7 


35.0 


34.0 


3.47 


3.32 


0.20 


0.30 


6.32 


6.32 


HD 210809 


09 lab 


21.2 


31.5 


32.0 


3.12 


3.31 


0.14 


0.20 


5.30 


5.80 


HD 30614 


09.5 la 


32.5 


29.0 


29.0 


2.99 


2.88 


0.10 


0.20 


6.04 


6.04 


HD 209975 


09.5 lb 


22.9 


32.0 


31.0 


3.22 


3.07 


0.10 


0.10 


2.15 


3.30 


HD37128 3 ' 


B0 la 


35.0 


28.5 
27.5 


29.0 i ' ) 


3.00 
2.95 


3.01 


0.10 
0.10 


0.10 


2.40 
3.01 


5.25 


HD 15558 


05 111(f) 


18.2 


41.0 


41.0 c) 


3.81 


3.81 


0.10 


0.08 


5.58 


7.10 


HD 190864 


06.5 III 


12.3 


37.0 


36.5 


3.57 


3.61 


0.15 


0.20 


1.39 


0.98 


HD 203064 


07.5 III 


15.7 


34.5 


34.5 


3.60 


3.60 


0.10 


0.20 


1.41 


2.58 


T Tl~"\ 1 At A n 

HD 191423 


09 III 


12.9 


32.5 


32.0 


3.60 


3.56 


0.20 


0.20 


< 0.41 


< 0.39 


HD46150 4 ' 


05 V((f)) 


13.1 


43.0 


40.0 


3.71 


3.71 


0.10 


0.10 


N/A 


1.38 


HD 15629 


05 V((f)) 


12.8 


40.5 


40.5 


3.71 


3.81 


0.08 


0.08 


1.28 


1.28 


HD5689 2) 


06 V 


7.7 


37.0 


36.0 


3.57 


3.66 


0.33 


0.20 


0.16 


0.17 


HD 217086 


07 Vn 


8.6 


36.0 


36.0 


3.72 


3.78 


0.15 


0.15 


<0.23 


<0.09 


HD 13268 


ON8 V 


10.3 


33.0 


33.0 


3.48 


3.48 


0.25 


0.25 


<0.26 


< 0.17 


HD 149757 


09 V 


8.9 


32.0 


33.5 


3.85 


3.85 


0.17 


0.17 


< 0.18 


< 0.15 


HD 149438 51 


B0.2 V 


5.3 


31.4 
31.9 


31.0 


4.24 
4.15 


4.00 


0.10 
0.12 


0.10 


0.009 
0.02... 0.06 


0.020 



Optical parameters tak en fr om 

'> Herrero et al]<2002l) 2) lHerrero et ail f2000) (unblanketed fastwind models) 

3) KudritzkietalJJl999l upper entries) and from Urbaneia (2004, lower entries, M scaled to R, = 35 R ) 

4) ^^^^^^^iT^^2t) (unblan keted plane-parallel H/He models) 

5) Kih^m^HuTfi99l|r and from lPrzvbilla & ButleJ l2004li with respect to wind properties (upper entries) 

and from R.M. (fastwind, lower entries); the limits of M correspond to velocity field exponents /? = 2.4. . .0.8. 
n) from a re-analysis by R.M. (fastwind) 
b) upper limit 

r) taken from optical analysis 



only IR spectra available we would favor this fit under the same 
circumstances anyway. 

Thus, in total we have four definite outliers (i.e., above the 
1-cr level and discarding HD 46150 because of its unclear sta- 
tus). HD 192639 and HD 209975 appear to be lower in gravity 
by -0. 15 dex when compared to the optical and the record hold- 
ers, Cyg OB2 #8C and HD 210809, give differences of -0.2 and 
+0.2 dex, respectively. In this sample, we might also include 
HD 46223 from Tab. |2] since the derived gravity is probably 
too low by a similar amount. 

Concerning mass-loss rates, the situation is as satisfactory 
as described above. First note that Fig. [2] also displays those 
stars for which we can only provide upper limits of M, and 
which we have compared. The only star missing in this com- 
parison is r Sco, however a comparison with both the "optical" 
mass-loss rate and the value cited in Tab.|4](which has been de- 



rived from an alternative IR analysis, cf. Sect. 16.3. 1> reveals a 
disagreement of a factor of roughly two (smaller and larger, re- 
spectively). Such a difference is not too bad, taken the intrinsic 
uncertainties at such low wind densities. We will come back to 
this problem later on. 

As expected from the non-linear increase of wind-emission 
as a function of M, the disagreement between optical and near 
IR mass-loss rates becomes smaller for larger wind-densities. 
For logM opt ^ -5.3, these differences are lower than 0.1 dex, 
which indicates the sensitivity of the Br r line wings (remem- 
ber that only the wings could be fitted in high M objects) and 
partly of Hen on this parameter. The differences obtained for 
the c orresponding equiv alent widths (observations vs. theory, 
cf. iLenorzer et al.ll2004l and Sect. 16.21 are thus almost exclu- 
sively due to the differences with respect to the line cores, 
which at present cannot be explained conclusively, though a 
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Fig. 13. Comparison of r e ff and log g derived from the optical and the near-IR. Asterisks, crosses and diamonds correspond 
to luminosity classes I, III and V, respectively. The displayed "error bars" correspond to maximum uncertainties claimed in 
optical analyses, namely +5% in T e g and +0.1 in logg. The two arrows in the left panel correspond to the objects HD 15570 
and HD 46150 (analyzed by means of unblanketed models in the optical) and indicate the average shift in positi on if blan keting 
effects would have been accounted for. The appropriate shifts in log g, applying the log g - T e g calibrations from Re polust et alJ 
(2004), have been indicated on the right, again by arrows, t Sco has not been included into the right figure, since the IR gravity 
could not be constrained due to missing Brl0/1 1. Concerning the remaining outliers, see Sect. [8] 
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Fig. 14. As Fig.[U] but for mass-loss rate M. Upper limits have 
been treated at their nominal values. The "error bars" corre- 
spond to mean uncertainties of +0.2 dex quoted for H a mea- 
surements, t Sco has not been included in this comparison (see 
text). 

relation to wind-clumping might be possible (cf. Sect.|9}. For 
the low M stars, on the other hand, no trend (particularly not 
towards considerably larger M) is visible. The above problem 
might be related to the physical conditions in the "intermedi- 
ate" wind, whereas close to the wind base no obvious differ- 
ence between the formation of H„ and Br y seems to be present. 

This hypothesis is consistent with the fact that the largest 
deviations are found at intermediate wind-strengths, i.e., in 
those cases where Br y is considerably refilled but does not 
show emission wings. The most prominent difference, A log M 
= 0.33, is found for HD 37128, if we compare with re - 
sults from the unblanketed analysis bv lKudritzki et alJ ( 1999). 



The corresponding re-analysis by means of blanketed models 
dUrbaneial 2004. see above) yields a moderately larger mass- 
loss rate 7 , and the remaining discrepancy, AlogM = 0.25 has 
been indicated in Fig.^J The difference for the other outlier, 
HD 203064, is of the same order, A log M = 0.26, and has to be 
considered as "real" as well, since all other parameters do agree 
(except for Yu e , which has only marginal influence on the de- 
rived mass-loss rate). For HD 217086, finally, differences occur 
only with respect to an upper limit, e.g., the IR analysis predicts 
a lower limit than the optical one. 

Concerning the helium abundance (not plotted), there are 
only two problematic cases. CygOB2 #7 has been discussed 
already in Sect. 16.3.31 and both the optical abundance (Yu e - 
0.2 ... 0.3) and the corresponding IR value (Yn e = 0.1 ) are 
uncertain, the latter due to the missing Hei lines (although we 
would derive a significantly lower T e g if the high abundance 
were true). The second outlier is HD 5689 with Yu e (optical) 
= 0.33 and FHe (IR) = 0.20. Again, this discrepancy is prob- 
ably irrelevant, since the optical analysis has been performed 
by means of unblanketed atmospheres, which are well-known 
to overest imate the helium abundance in a number of cases (cf. 
Repo lust et all2 004, particularly Sect. 7.2). 

8.1. Comments on hydrogen collisional cross sections 

In Sect. I3.ll we briefly discussed the importance of consistent 
collisional data for the resulting IR li ne profiles. We outl ined 
recent calculations performed by Przv billa & Butlerl 12004) and 
references therein. These authors provide a recipe for an "op- 
timum" choice of collisional data, based on a number of com- 

7 though at a lower effective temperature and a somewhat larger 
value of /?, which, in combination, are canceling each other with re- 
spect to M. 
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Fig. 15. Influ ence of different collisional cross sections for hydrogen: some prototypical examples. Solid: Best fit with cross 
sections from Giovanardi et al. ( 1987), as in Figs.|9]to^2 Grey: Models with identical parameters, but cross sections calculated 
according to Przvbilla & Butler (2004). See text. 



parisons with observations, comprising mostly BA-type dwarfs 
(including t Sco) and supergiants, whereas only one test has 
been presented for an O-star, the 03.5 dwarf HD 93250. 

After incorporating their data into our version of fastwind 8 , 
we have subsequently tried to analyze our observed dataset. 
First let us mention that this modification gave rise to changes 
in the hydrogen lines alone, since in all cases the tempera- 
ture structure (being dependent on the collisional bound-bound 
rates) remained almost unaffected, with maximum changes on 
the order of 100 Kelvin. 

Unfortunately, however, it turned also out that, again in all 
cases, the hydrogen line cores became stronger (in agreement 
with the findings by Przybilla & Butler), whereas the line wings 
are barely affected, as shown in Fig. [EJfor some prototypical 
examples taken from the fits in Figs.|9]toE] In a few cases this 
might actually lead to an improvement of the situation, e.g., 
for HD 46223 or Cyg OB2 #8C, which actually need gravities 
higher than those derived above. (Remember that increasing 
the gravity results in shallower line cores, cf. Sect. 14. 1> . Since, 



8 Note that Przybilla & Butler have used a similar version to check 
part of their calculations. 



however, the IR-graviti es bas ed on our standard collisional data 
from Giovanardi et al 1 lll987l) were found to be consistent with 
the optical ones in the majority of cases, models based on the 
alternative data by Przybilla & Butler would consequently lead 
to an overestimate of gravities. 

The same would be true for the mass-loss rates. Using the 
new data would sometimes improve the situation, e.g., any 
central emission inside Br y (if present) becomes reduced, cf. 
HD 46233 and t Sco in Fig. ^] Actually, this reduction is 
the origin of the lower m ass-loss rate of r Sco as derived by 
Przvbilla & Butler (2004): If M is decreased, the strength of 
this feature increases again, at least if the winds are very weak, 
due to subtle NLTE effects and unrelated to any direct wind 
emission. 

For "normal" winds, on the other hand, where wind emis- 
sion plays the primary role, the deeper cores predicted by the 
"new" models would necessitate higher mass-loss rates. Thus, 
the present situation would get worse, at least in those cases 
where M is no longer derived exclusively from the wings. As 
an example, consider HD 203064 in Fig. which presently 
has an IR mass-loss rate which is a factor of two larger already. 
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Using the new dataset would further deteriorate this discrep- 
ancy. 

Insofar, the preference for our standard set of collisional 
data is triggered solely from the results described above, 
namely from the generally satisfactory agreement between the 
IR and optical analyses/or those objects analyzed in the present 
investigation. We do not argue that one set or the other is bet- 
ter, but point out that in these cases our standard dataset gives 
results which are more consistent with the optical. 

Of course, we have also looked into some of the details 
responsible for the differences obtained. It turned out that the 
NLTE departure coefficients are astonishingly similar, when 
comparing the results from both collisional datasets. There are 
only (very) subtle differences in those regions where the line 
cores are formed giving rise to the deeper profiles if the dataset 
by Przybilla & Butler is used. Either the lower level of the tran- 
sition is slightly more populated, or the upper one is slightly 
less populated. The obvious discrepancies are then induced by 
the extreme sensitivity of the IR line formation on such subtle 
differences. As has already been argued about the formation of 
Hei lines (Sect. I4.2L this discrepancy would barely be visible 
if the lines were situated in the optical. Insofar, not only the 
data but also the numerical treatment plays a crucial role. In 
any case it is quite astonishing how well the observed profiles 
can be simulated. 

9. Summary and conclusions 

In this paper, we have analyzed 25 Galactic O and early B-stars 
by means of H and K band spectroscopy. The primary goal of 
this investigation was to check to what extent a lone near-IR 
spectroscopy is able to recover stellar and wind parameters de- 
rived in the optical. This is critical to our desire to precisely an- 
alyze the hot, massive stars, deep within the disk of our Galaxy, 
and in particular the very young, massive stars just emanating 
from their birth places. 

Most of the spectra have been taken with subaru-ircs, at an 
intermediate resolution of 12,000. In order to synthesize the 
strategic H/He lines present in the H/K band, we have used 
our recent, line-blanketed version of fastwind. In total, seven 
lines have been investigated, three from hydrogen, including 
Br r serving as a diagnostic tool to derive wind-densities, two 
Hei and two Hen lines. For two stars, we could make additional 
use of Hei2.05 (singlet) which has been observed with irtf- 
cshell. Apart from Br y and Hen2.18, the other lines are pre- 
dominately formed in the stellar photosphere, and thus remain 
fairly uncontaminated from more complex physical processes, 
particularly clumping. 

In our attempt to prepare all required broadening functions, 
it turned out that at present we have to rely on the Griem ap- 
proximation for Stark broadening (important for hydrogen and 
Hen), since the corresponding published data (based on the 
more exact VCS approach) suffer from numerical problems, 
particularly for the members of higher series. 

First we investigated the predicted behaviour of the strate- 
gic lin es, by means of a large model grid described in lPuls et"afl 
(2005). Interestingly and in contradiction to what one expects 
from the optical, almost all photospheric lines in the H and K 



band (from H, Hei and Hen ) become stronger if the gravity de- 
creases. In Sect.0] we have carefully investigated the origin of 
this rather unexpected behaviour. 

Concerning H and Hen, it is related to the particular be- 
haviour of Stark broadening as a function of electron density, 
which in the line cores is somewhat different for members of 
lower and higher series. For the latter, the cores become deeper 
when the density decreases, and contribute more to the total 
line strength than in the optical. 

Regarding Hei, on the other hand, the predicted behaviour 
is due to some subtle NLTE effects resulting in a stronger over- 
population of the lower level when the gravity decreases, so 
that the source function becomes weaker and the profile deeper, 
i.e., stronger. This strong dependence of the profile on the 
source function is a direct consequence of the IR line formation 
with hv/kT « 1. If those lines were situated in the optical, on 
the other hand, optical depth effects would dominate, leading 
to a decrease of line strength due to a lower number of Hei ab- 
sorbers. This explains the different (and "normal") behaviour 
of, e.g.,Hei4471. 

In Sect. we ha ve compared our calcu lations with results 
presented recently bv lLenorzer et alJ (|2004j), utilizing the alter- 
native NLTE model atmosphere code cmfgen. In most cases, 
we found reasonable and partly perfect agreement. Only the 
Hei2.05 singlet for mid O-types suffers from some discrepancy, 
in ag reement with our analogous findings for optical Hei sin- 
glets dPuls et al.l2005l) . 

After carrying out the analysis for our sample described 
above (and in agreement with the predictions from our model 
grid), we find that an H/K band analysis is able to derive con- 
straints on the same set of stellar and wind parameters as it is 
known from the optical, e.g., T e fi, \ogg, Yn e and optical depth 
invariant Q, where the latter yields the mass-loss rate M if stel- 
lar radius and terminal velocity are known. For cooler objects, 
when Hen is missing, a similar analysis might be possible if 
Hei2.05 (provided to be described correctly) is available and 
the helium content can be adopted (due to the almost orthog- 
onal reaction of Hei2.05 and Hei2.11 on r e ff and logg). This 
should be possible for very young objects containing unpro- 
cessed material. 

For future purposes, when no UV observations will be 
available, the terminal velocity Voo has to be taken from cali- 
brations, as it is true for the velocity field exponent (3, at least 
in those cases when no emission wings in Br r are visible. 
Concerning the determination of R„, a similar strategy as in 
the optical might be developed, utilizing infrared colors and 
distances. 

For most of our objects, we obtained good fits, except for 
the line cores of Br r in early O-stars with significant mass-loss 
(see below), and except for the fact that particularly at mid O- 
types Br 10/ 11 could not be fitted in parallel. We have argued 
that this discrepancy is similar t o the problem in th e optical, 
concerning Hen4200/4541 (e.g.. iHerrero et alJEoO^ . Due to 
the similarity in the involved levels and broadening functions, 
we have speculated about a possible defect of these broaden- 
ing functions for transitions between members of higher se- 
ries. The largest discrepancy, however, was found for the line 
cores of Br y . First note that this problem is not particularly re- 
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lated to our code, since also cmfgen exhibits the same short- 
coming. Whereas the observations show Br y mostly as rather 
symmetric emission lines, the models predict a P Cygni type 
line, with a comparably deep core which is never observed. 
Note that this type of profile can only be created if the ratio of 
departure coefficients for the involved levels (n = 4 — >7) devi- 
ates strongly from unity (cf. IPuls et al"1l 19961) . whereas a ratio 
close to unity would just give the observed symmetric emis- 
sion profile. One might speculate that this can be achieved due 
to a stronger influence of collisional bound-bound processes, 
which, e.g., might be possible in a strongly clumped medium. 
Remembe r, however, that Br v typically forms inside the Hq,- 
"sphere" jLenorzer et al.ll2004T) . where the de gree of clump- 
ing is usually thought to be moderate (e.g., Ma rkova et alJ 



2004; Repolust et al. 2004), though a recent investigation by 
Bouret et alJ (| 2005) strongly indicates the opposite. 



After having derived the stellar and wind parameters from 
the IR, we have compared them to results from previous opti- 
cal analyses, in an almost strictly differential way, since most of 
these results have been obtained also on the basis of fastwind. 
Overall, the IR results coincide in most cases with the opti- 
cal ones within the typical errors usually quoted for the corre- 
sponding parameters, i.e, an uncertainty in of 5%, in logg 
of 0. 1 dex and in M of 0.2 dex, with lower errors at higher wind 
densities. In most of the cases where we have found discrepan- 
cies beyond these errors, their origin could be easily identified. 
Definite outliers above the 1-cr level where found in four cases 
with respect to log g and in two cases for M, at intermediate 
wind-strengths. Given the 1 -<x character of the discussed errors, 
these mismatches are still not too worrisome, since all identi- 
fied discrepancies lie well below 2-<x, and no trends are visible. 
Particularly with respect to logg it is quite possible that er- 
rors exceeding the nominal error of 0. 1 dex do exist in certain 
cases, both in the optical and the IR, relat ed to fits whic h are 
not at their (global) optimum (e.g., Moki em et alJ Eo05). Re- 
assuring is the fact that the mean difference between NIR and 
optical gravities is very close to zero, Alogg = -0.025 + 0.1. 

As a by-product of our investigation, we could determine 
the (IR-) stellar parameters and the log Q value for three 
dwarfs, which have not been analyzed in the optical so far. 
Though we have shown that the derived quantities are rather 
insensitive to uncertainties in R, and Voo, these objects need to 
be checked in the optical. 

Let us highlight one additional "bonus" obtained from the 
infrared. In those cases when a star has an extremely weak 
wind and the core of Br r can be resolved (requiring a very 
low rotational speed), the central emission will give us a clue 
about the actual mass-loss rate and not only an upper limit, 
as is true for the optical. An example of this kind of diagnos- 
tics is t Sco. Particularly with respect to recent investigations 
of yo ung dwarfs with surprisingly weak winds (M artins et alJ 
120041 this will turn out as an invaluable source of informa- 
tion (even more, if coupled with observat ions of Br a , e.g. 
iNaiarro et all 1 9981: IPrzvbilla & ButleJ2004T) . 

After finishing this investigation, we are now able to con- 
strain the observational requirements to perform such a de- 
tailed IR-analysis (but see below). Most important is a high 
S/N, because most of the lines to be investigated are extremely 



shallow, at least for the hotter objects. The resolution must be 
sufficient to disentangle the line cores from the wings (partic- 
ularly important for Br y ) and to obtain reasonable clues about 
any contamination due to reduction problems. As for the re- 
quired set of lines, almost all lines analyzed in the present pa- 
per are necessary to obtain useful constraints, maybe except 
for BrlO, since Brl 1 seems to be less contaminated. Since both 
Hen lines behave very similarly and show the same degree of 
consistency or disagreement (if present), one of those two lines 
might be discarded as well. 

In the last section of this paper, we have argued that our 
standard implementation of hydrogen collisional cross sections 
seems to give results which are in better agreement with the 
optical results for our sam ple of hot objects, compar ed to the 
data suggested recently bv lPrzvbilla & Butlerl J2004I) . In view 
of their findings, namely that for cooler stars their prescription 
gives more consistent results, this discrepancy has to be clar- 
ified in future work. The same of course is true regarding the 
severe mismatch of the Br r cores. A first step will require to 
include clumping and to investigate to what extent this process 
might improve the situation. 

The value of a reliable quantitative analysis for hot, mas- 
sive stars based entirely in the infrared cannot be overstated. 
Most obvious, it will allow the evaluation of massive star char- 
acteristics at an evolutionary stage significantly earlier than has 
ever been possible before. The influence of disk emission may 
render the photosphere of some very young massive stars in- 
accessible. We suspect, however, that among the most massive 
stars, around mid-O or hotter, the disk will be destroyed well 
before even near-infrare d studies would be feas ible due to the 
very short disk lifetime (Wats on"& Hanson! 19971) . 

A combined H and K band analysis with the S/N we present 
here won't always be possible, nor is it required. A S/N of just 
150 and resolution of order R ~ 5000 will typically be suffi- 
cient for a reasonable an alysis, provide d the star is not a slow 
rotator (see discussion in lHanson et al.lEo05l) . How deep into 
the Galaxy one can probe is a strong function of the line of sight 
extinction and the intrinsic stellar brightness. For late-O dwarfs 
with Ay ~ 20, one can only probe a few kiloparsecs. However, 
for early-O dwarfs or most O and early-B supergiants, even 
with Ay ~ 20, one should be able to probe the entire solar cir- 
cle. For this level of extinction, it is the //band measurements 
which pose the greatest challenge, though S/N ~ 120 - 150 
should still be attainable with an 8-m class telescope. For still 
redder stars, the H band would be too dim for useful spectra. 
Once m.K > 14-15 (OB supergiants at the far side of the Galaxy 
with Ay ~ 30), the corresponding low S/N for such a dim star 
wouldn't allow a robust quantitative analysis, though spectral 
classification is still possible. Thus, excluding sight-lines where 
the extinction has become extreme (Ay > 30), nearly every lu- 
minous OB star within our Galaxy becomes accessible through 
near-infrared spectroscopic studies. 

One of the powers of quantitative analysis is its ability to 
determine absolute magnitudes. When OB stars are in clus- 
ters, those cluster distances will be robust, giving us clues to 
the structure and nature of the presently, poorly understood, in- 
ner Milky Way. OB stars serve as secondaries to massive com- 
pact objects. Because the extinction is typically high for such 
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systems found in the inner Galaxy, a NIR analysis of the OB 
companion provides the only means for making critical mea- 
surements to constrain these fascinating systems. In truth, be- 
cause many more O stars within our Galaxy are visible in the 
NIR than the optical by almost two orders of magnitude, the de- 
velopment of a robust quantitative analysis in the infrared will 
stimulate entirely new, important results on massive stars, their 
formation and evolution and numerous valuable insights into 
the inner workings of our Milky Way Galaxy. 
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